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GASDYNAM1C  LASERS:  THEORY, 

EXPERIMENT  AND  THE  STATE-OF-THE-ART 
by 

John  D.  Anderson,  Jr. 

Professor  and  Chairman 
University  of  Maryland 
College  Part,  Maryland,  U.S.A.  20742 

I.  GAS DYNAMIC  LASERS:  THEORY  and  EXPERIMENT 

A.  introduction 

Consider  in  your  inagination  a supersonic  wind  tunnel,  such 
as  those  cornonly  found  in  many  aerodynamic  laboratories.  However, 
when  this  particular  tunnel  is  turned  on,  we  do  not  measure  the  lift 
and  drag  on  an  aerodynamic  model,  o.  the  pressure  and  heat  transfer 
distributions  over  a surface.  Instead,  when  the  switch  is  thrown 

for  this  particular  tunnel,  we  see  a^rijiower^ 

from  the  test  section.  Indeed,  this  is  not  a wind  tunnel  at  all.  but 

rather  it  is  a gasdynanic  laser. 


Fig.  1 


Schematic  Drawing  c>f  G.is  dynamic 
Laser 


As  shown  sc"  eniatically  in  Figure  1,  a gasdynamic  laser  takes 
a hot,  high  pressure  mixture  of  gases  (usually  C02,  N2  and  H20  or  He) 
and  expands  this  mixture  very  rapidly  through  a supersonic  nozzle.  During 
the  expansion,  the  gas  is  turned  into  a laser  medium  (a  population 
inversion  is  created).  The  supersonic  laser  gas  then  passes  into  the 
test  section  (laser  cavity),  where,  if  mirrors  are  placed  on  both 
sides  of  the  test  section,  a bean  of  laser  energy  is  extracted  perpendicular 
to  the  flow.  The  supersonic  stream  then  enters  a diffuser,  where  it  is 
shocked  down  to  subsonic  speeds  and  generally  exhausted  to  the  atmosphere. 

Why  are  gasdynamic  lasers  so  important?  The  answer,  purely 

and  simply,  is  that  they  are  very  hi  eh  power  devices.  In  fact,  gasdynamic 

lasers  (GDL's)  have  spearheaded  a breakthrough  in  high  energy  laser  devices. 

For  example,  a multimode,  continuous  wave  power  output  of  63  KW  from  a 

C02-,’J2  GDL  has  been  reported  by  Gerry1,  and  an  average  laser  power  of 

400  KW  has  been  extracted  for  four  milliseconds  from  a shock  tube  GJL 

2 

by  Klosterman  and  Hoffman  working  in  the  laser  group  at  tiie  University 
of  Washington.  A major  factor  in  this  breakthrough  is  that,  unlike 
early  electric  discharge  lasers  with  their  attendent  problems  of  arc 
discharges  in  large  volumes,  GDL's  can  be  sealed  to  large  sizes  without 
major  physical  complications.  A picture  of  a large-scale  GJL  is  shown 
in  Figure  2. 
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Fig.  2 Gtlc/y^c  /»«--  ^ 

A historical  note  is  in  order  hero.  The  concept  of  a 
asdyn.r  TcT^was  initiated  as  early  as  1962  by  Basov  and  Oraevstii3, 

10  suggested  tint  population  inversions  in  molecular  systems  could  be 
re.  ted  by  rapid  heating  or  cooling  of  the  system.  Subsequently,  in 
900  Hurl e and  ilertzberg4  suggested  that  such  cooling  and  popu’ation 
n versions  could  be  obtained  in  the  rapid,  nonequilibrium  expansion  of 
r.  initially  hot  gas  through  a supersonic  nettle.  They  considered  the 
.pacific  case  of  electronic  level  inversions  in  expanding  So.  but  were 
msuccessful  in  measuring  suen  inversions  in  the  laboratory.  Then, 
tantowitz  combined  this  idea  with  his  previous  work  on  vibrational 
lonequilibrium  in  C02  (see  Reference  5),  and  in  1966,  he  and  a group 
of  physicists  and  engineers  at  the  AVCO  Everett  Research  Laboratory 
operated  the  first  GUI,  using  a mixture  of  (sen  Reference  1). 

At  about  the  same  time.  Basov  ct  al6  carried  out  a theoretical  analysis 
of  population  inversions  in  C02-!l2  expanding  mixtures,  and  predicted 
that  substantial  Inversions  can  indeed  occur  under  suitable  conditions 
for  the  90S  mixture  r.-.io  and  the  nozzle  reservoir  pressure  and  tempo, turn. 
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However,  ,,asov  missed  an  essential  point,  namely,  that  the  inclusion 


of  a "catalyst"  such  as  H2O  or  He  is  necessary  for  the  attainment  of 

reasonable  population  inversions.  Easov  used  older  values  of  the 

vibrational  r-ergy  exchange  rates  winch  resulted  in  optimistic  values 

of  the  population  inversions;  indeed,  using  more  recent  rates  as  compiled 

7 8 

by  Taylor  and  Bitterman,  Anderson  har  shown  that  the  mixtures 

considered  by  Basov  et  al  do  not  yield  significant  inversions.  The 

group  at  AVCO  recognized  this  fact,  and  have  reported  experiments  using 

CCk-llg-lUO  mixtures  as  early  as  1966^.  Theoretical  calculations  showing 

8 

the  beneficial  role  of  I^O  were  reported  by  Anderson  n 1969.  The  exact 

role  of  ILO  or  lie  as  a catalyst  is  discussed  later  in  these  notes. 

Since  these  initial  experiments  and  calculations,  the  technology 

and  fundamental  understanding  of  GuL's  lias  grown  precociously. 

9-12 

This  growth  is  exemplified  by  experiments  carried  out  in  arc  tunnels  , 

shock  tunnels2’^’^"^  and  combustion  driven  devices 1 ’ ' 

and  by  theoretical  calculations  reported  in  References  6,3,11,12,21-26. 

Moreover,  gasdynanic  lasers  using  CO  as  the  lasing  medium  have  also 

been  demonstrated  . It  is  the  purpose  of  these  notes  to  bring  the 

reader  up-to-date  in  this  state-of-the-art. 


B.  ELEMENTARY  PHYSICS 

1 . Enemy  Levels  and  Population  Inversion 

Consider  a collection  of  molecules  in  a gas.  Pick  one  of  the 
molecules  and  examine  it.  The  molecule  is  moving  about  in  space  — hence 
it  has  translational  energy;  it  is  rotating  about  its  principal  axes  — 
hence  it  has  rotational  kinetic  energy;  the  atoms  that  make  up  tne  molecule 
may  be  vibrating  back  and  forth  from  some  equilibrium  position  --  hence 
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It  has  vibrational  kinetic  and  potential  energy;  and  finally,  the 
I electrons  irove  abet  the  nuclei  of  the  molecule  - hence  it  has 

' kinetic  and  potential  energy  of  electtonic  origin.  However,  the 

I amazing  quality  of  these  various  form  of  molecular  energy  is  that 

| they  can  not  be  any  arbitrary  value.  Rather,  one  of  the  most 

important  discoveries  of  modern  physics  is  that  the  molecule,  at  any 
; given  instant,  has  to  occupy. . 

* levels,  i.e.,  the  energy  values  of  a molecule  are  chastized.  This 

I is  shown  schematically  in  Figure  3.  Mow,  instead  of  one  molecule, 

imagine  the  whole  collection,  say  106  molecules,  and  look  at  the  first 
1 quantized  eh.  gy  level  c„  (the  ground  state).  Thera  may  be  400.003  molecules 

in  this  level  at  some  instant  in  time.  The  number  430,000  is  called  the 
population  :io  of  the  c0  level,  ilext,  look  at  the  first  exited  level. 
t } there  may  be  200,000  molecules  in  this  level,  hence  II,  « 200,000  * 
the  E:  eolation  of  the  c,  level;  and  so  forth.  The  set  of  numbers, 

H , tl2,  ...  H,,  ...  It  called  the  gH»1rt<»»  dlstrihutla  »»«•  tne 
ene  igy  levels  of  the  gas.  The  nature  of  this  population  distribution 
is  of  vital  importance  for  loser  action.  For  example,  consid.r 
vibrational  enemy  of  a molecular  gas.  If  the  gas  is  in  thermodynamic 
equilibrium,  the  population  distribution  will  exponentially  decrease 
with  incr  .sing  c,.  as  shown  on  the  left  of  Figure  4,  i.e..  it  will 
follow  a r tzrm.nn  distribution,  ft  major  characteristic  of  this  equil  ibria.' 
distribution  is  that  IS,*,  < II,.  H°"'=vor,  if  the  gas  is  disturbed  at 
any  instant,  say  by  moans  of  an  electric  discharge,  or  by  a very  sudden 
temperature  change,  than  the  population  distribution  can  baco.ua  a 


noncguilibrium  distribution,  .nd  it  is  even  possible  that,  at  some  instant. 
„ This  situation,  where  the  number  of  molecules  in  a hinher 

1 y i ri energy  level  is  orator  than  the  number  in  a lower  lyjsa  u'vel  * 
is  called  a ^uUti^in^AS!!-  population  inversion  is  the  essence 

of  laser  action,  and  the  attainment  of  this  population  inversion  by  means 
of  rapid  cooling  of  the  gas  is  the  essence  of  a gasdynaiiic  laser. 
Specifically,  when  a gas  expands  through  a supersonic  nettle,  as  shown  . 
in  Figure  5.  the  gas  temperature  decreases  very  rapidly.  Indeed,  a fluid 
element  moving  through  the  nontie  c,  experience  temperature  changes 
as  high  as  IS6  degrees  per  second  in  the  throat  region  of  the  nettle, 
in  fact,  expansion  through  a minimum-length  supersonic  nozzle  is  tie 
fastest  practical  me  ns  of  cooling  a gas.  Therefore,  it  is  reasonable 
to  expect  that,  for  certain  conditions,  a population  inversion  can  be 

created  in  such  ar  expansion. 


Fig.  5 
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Witi)  the  above*  discussion  'iri  mind,  tin:  pivotal  and  first 
problem  in  studying,  understanding  and  analysing  a gasdynamic  laser 
is  to  examine  the  nonofrilibrim  flow  tiv.  nugh  a nozzle,  and  the 
behavior  of  this  flow  as  it  moves  downstream  through  the  laser 
Cdvi\_, , opec i f i cal  ly . wo  wish  to  calculate  the  population  distribution 
as  it  varies  through  the  nozzle  for  a mixture  of  C02-ri2  and  11^0  or  He, 
and  we  wish  to  examine  the  conditions  under  which  a population  inversion 
will  develop.  This  will  be  the  subject  of  Section  C of  these  notes. 

However,  before  we  address  this  problem,  we  should  answer  the  .allowing 
question . 

2.  f'hy  Poes  a Po:v|  i at  ion  ’rivers  ion  Make  a Laser  Work? 

Cons  ider  two  energy  lcr.  jls  of  the  molecules  in  a gas.  There 
are  throe  mechanism:;  by  which  there  energy  levels  interact  with 
radiation:  spontaneous  emission,  stimulated  emission,  and  absorption. 

These  are  illustrated  in  Figure  C. 

Now,  consider  a single  molecule  in  a field  of  radiation, 
whore  the  radiative  intensity1  per  unit  frequency  is  I 

v" 

( * *2 1 ) v - nis'i.e.r  of  spontaneous  red  stive  trar.itions  p.-r  second  per 
particle  per  unit  frequency . This  is  called  the  Einstein 
spun  tar,?  ous  emission  coefficient,  and  sometimes  the  Einstein 
t r a n s i t i o n p r o b a t > i 1 i ty . 

(321)vIv.  number  of  induced  radiative  transitions  per  second  per  particle 

per  unit  fro  money.  (B^  )v  is  the  Einstein  induced  "emission  coefficii 

(^1 2^ v^v'  11  ,or*  Cji%  >*>■’- di alive  Irons i ti onr.  per  second  per  particle  per  unit 
frequency  due  to  absorption.  (3, »,)  is  the  Einstein  absorption 
coc  f i i ci  oil  i . tc  v 

Consider  a slid)  of  gas  of  geometric  thickness  dx,  as  shown 
in  Figure  6.  Let;  1 !,o  the  radiative  intensity  (energy  or  second  per 

uni  area)  incident  on  the  left  side  of  the  "Tub,  and  let  ( c!I ) bo  the 

BEST  AVAILABLE  COPY 
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change  in  intensity  while  traversing  the  slab.  Also,  let  and 
be  the  nu:;ber  of  molecules  per  unit  vole -e  in  the  and  c ^ energy 
levels  respectively.  Then,  from  t-e  above  definitions,  we  can  simply 
write  down 


( c/Z - [{ M -/•  ( 3*/  ),~  M Ml  { ^2  )r.'Z\rMj'  ‘ >~*M  ( * ) 


i/V.i  fc/ncuS 
t'rru  si'/i’/i 


/>i ,/.  . t . / 


C)bsc.-y  //  C /-> 

For  the  infared  wavelengths  characteristic  of  mol  ecu’ ..r  lasers,  the 

12  9] 

spontaneous  emission  is  negligible  . Also,  in  general,  (B^)^  g~*  ‘^I2^v* 
as  shown  in  Reference  12  and  in  most  modern  physios  textbooks.  Here, 
g^  and  are  statistical  weights  f the  levels.  Hence,  eg.  (1) 

l>  omer. 


(c/T\ 


2~ 


(%*  (Jf,  /%-*/)  J‘l'~  c/* 
/ 


(2) 


How  we  ccn  answer  our  original  question.  Looking  at  eg.  (2),  (dl)  / 1 

_ V V 

9i 

will  be  positive  when  (— Ho-iU>D,  i.c.,  when  a popular'  'n  inversion 

«2  c 1 

exists  between  the  r.n  and  c-|  energy  levels,  hence,  a population  inversion 
1 ea  ns  to  an  ar  • ”!  i fic  a tion  of  the  radl : Live  inter  si  tv  I . This  is  the 
laser  off  ct.  and  this  is  why  a population  inversion  rakes  a laser  work.! 


3.  Definition  and  Cal  elat  ion  of  S:  all -Signal  Gain 


The  small-signal  laser  gain  is  an  extremely  important  figure 
of  merit  for  gas  las.r  devices;  it  is  to  laser  physicists  what  lift  and 
drag  coefficients  are  to  aerodynamicists.  It  is  also  a direct  measurement 
of  the  papulation  inversion.  Ev  ry tiling  else  being  equal , the  higher  trie 
small-signal  gain,  the  easier  leer  action  can  be  obtained  in  a gas. 
Therefore,  because  we  h.ve  just  discussed  the  reasons  why  a population 
inversion  k -s  • laser  wor!: , we  now  extend  this  discussion  to  the 
analysis  of  snail -sign  1 ;;jin.  Also,  th  following  results  are  necessary 


11 
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for  our  subsequent  discussions  on  ihe  theory  end  experimental 
measurements  of  gasdynan'c  laser  performance. 


i.et  us  consider  stimulate.-'  emission  and  absorption,  ignoring 
spontan  ous  emission  as  stated  previously,  Consider  again  the  slab  of 


gas  of  geometric  thickness,  dx, 

{:/!)  t-/fvE 

with  radiation  of  intensity  I 

A (per  unit  frequency)  iucid  nt 

normally  on  the  slab  as  shown  in 
Fig.  6.  The  ri  Motive  intensity 
absorbed  in  a nivan  rnsctr  1 line 

c 1 y— 

of  the  gas  is  (dl)v  (per  unit  frequency) 
Due  to  line  broadening  effects  in  the  gas  (natural,  Doppler  and  collision), 
(dl)  varies  over  a narrow  frequency  rauoe,  as  shown  in  the  sketch  on 

V 

the  left.  This  one.  rise  to  toe  1 i r: :•  sha;  - The  intensity  of  radiation 
between  v and  v+dv  is  (dl)  dv.  The  integrated  radio  vc-  intensity 
(total  intensity  absorbed  by  t:  line)  is 

d. 

The  spectral  absorption  coefficient,  , is  define  I 

(dl)  = -a  I dx  (3) 

v 7v  v v 

Examining  cq.  (3),  if  I d.:crear-'  as  it  traverses  the  slab  of  gas,  then 
dl  is  negative  and  a is  a positive  number.  On  the  other  iund,  if  I 
increase"  as  in  the  laser  effect,  then  dl  is  positive,  and  &v  has  to  be  a 
negative  quantity.  Rather  than  cl;  a 1 with  a "neg  Live  rbsoi  ,.-t  i on  coe  > icier.  L' 
a nc',;  coefficient  is  defined  os 


(dl)  = G I dx 

' ' V 0 V 

V 
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where  Gq  is  the  small -signal  1 in  coefficient. 


V“ 


(5) 


Since  the  snail  signal  g.  a coefficient  is  the  same  as  a "negative 
absorption  coefficient",  let  us  pursue  orr  discussion  in  terms  of  the 
absorption  coefficient.  Comparing  c.  ;s.  (2)  and  (3),  we  sec  that 


= hr  (Q*\  (a/,  - ^ /14  ) 


(0) 


It  is  com  an  to  express  (B.J  in  toms  of  the  Einstein  coefficient 

It  v 

for  spont  -neous  enis  ion,  ( 1 ) v * wiv'rc  (seo»  for  example.  Ref.  „3) 

n / 1 

S n / "r~ 


( &:/)y  = 


( )v 


(*/) 


Thus,  e ...  (G)  beccr.es, 

Eq.  (7)  gives  the  spectr  1 absorption  coefficient,  a . Define  the 
into  nntrd  absorption  coefficient  as  ( ct^dv  where  the  integral  is  taken 
over  the  entire  spectral  line. 


Thus 


= £}(. A;-  §/:•;)/ V*4 


/-  - -v 

Rote:  The  interval  of  int  gr-  ion  over  v is  snail  (lines  are  usually 

0 9 

narrow),  so  that  v in  the  factor  c /oW  in  eq.  (ft)  can  be  con: icie>v  .cl 
constant.  Letting  j" (;'21 ) vrfv  ^ ^2i  ’ i;c  heVl-‘ 


(2) 


y-  ( /<  " 


C ' 

&7\' 


:,/?v  (*!-  §'  M ) 


C) 
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Sin  ; - = X,  cq.  (9)  becomes 

- -& /;-  ^ 1 ; <iu> 

V 

Eq.  UO)  gives  the  it;t  ir'U'1  absorption  cocffic.  .out.  Eq.  (10)  is  nore 
useful  then  cq.  (7)  because  , n_o_t  (A?1)v,  is  available  from  laboratory 
measurements.  To  sea  this  mare  clearly,  recall  eq.  (2). 

(cir)r  = -fax.  (a;- 

( lYJ)  v = ^ ~J  ( SU/)s-  ( A't  — /e  )-~r  A >-  c/s 

Ur)>-  " £fa'!»)rfa-p4)Z-c'y  m 

How,  i nt.e,; rate  cq.  (11)  over  t tiii^ncss  of  the  slab.  L t L - slab 
tlvichn.  s /'t'l.f  -j  „ j 


c i.  r 

ru  _ 


> ‘ / •'  ' (/ 


K'nf-  ....  .sVr  ^ ~ - fa  ) L 

— - “ C ' 

Jr 

*w 

For  small  values  of  the  exponent,  i.e.,  sie.ll  L, 

P-  = /-  ■' — ( 

in  (A  ' 

The  total  radiative  intensity  -bsorhcrl  over  the  entire  spectral  line  is 


•abs*  v;l,ore 


j'-,,  = / ( A'.  - Tv- , ) r/\ 

On/  A 


Assume  that  a nt  i r.  1 ! ,t  : *ureo  is  us  •’  in  an  . sorption  experv  nt, 


I «i  i c fi  , ■ }•  r i ■ ■ . 


r i f v ii  aval  of  tit*-  s,  ' st  eal  I it 


Confining  cqs.  (Ik)  and  (13): 

X/,  * /_  Lj ~x,„  +Tr.,„  ('h,K  (V~§1  "*)£*'*- 


u/.j-  - r • 


/// 


Tjt"  ^ /y'  ~ ~fxA>')  L r//^ 


y~  ...  2l7:Z--/  / / //  \ / 7-  , 

-4^-  iv  <-7lv  pAA/xJ*~ 


I.l.  is  measured  in  the  lab.  In  turn,  A.,,  is  calculated  from 

eq.  (14).  Therefore,  it  is  U] , not  (A^)v>  that  is  usually  .obtained 

from  experiment.  Thus,  os  stated  earl  ;or,  eq.  (10)  Ur  the  integrated 

absorption  coefficient  is  more  directly  related  to  experimental  observe*. ion 

than  is  eq.  (7)  for  the  spectr.';!  absorption  coefficient.  Also,  keep  in 

mind  that  A^  is  a physical  constant  of  the  molecule. 

We  have  seen  that  in  a conventional  absorption  experiment 

using  a co:  tiircus  light  source,  it  is  the  integrated  absorption 

coefficient  f n riv  and  hence  A-,,  which  is  obtained  from  the  dat  . 

Values  for  A^jt  and  also  for  the  radiative  lifetime, 

. 1_ 
hi 

arc  comr.orily  quoted  in  the  1 itorature—  

Now,  consider  l_as.  rad iation.  This  radiation  is  highly 

monochromatic.  The  laser  cavity  is  generally  tuned  to  emit  radiation 

over  a very  small  frequr  y interval.  This  interval  is  usually  much 

smaller  than  the  liuewidlu  of  a spectral  line, 
of, 


Li  lie  !.//«-.  f>t-  u /'  r\ 
/’I- a i * a/ 1'  ii  i*  </  b y T 11 
////  <- 


k Ou  ll  ^ £ ".li  itoV  fa. 
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Li 
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0 
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J 

t . 
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the  spectral  absorption  coefficient,'^. 


Nov/,  lot  us  relate  to  the  integrated  absorption  coefficient 


Define  a line  sh  .pe  factor, 
g(v,vo)  as: 

/"<-  rU-  05) 


where 


is 


obtained 


the  interval  of  integration  is  over  the  entire  spectral  line. 

For  pres-  ares  above  10-20  tor r, the  line  shape  is  dominated 

by  coll  i r- i c;i  broadening  (or  nr  * sure  1 ~o'  do.vinn) . It  has  been  demonstrated 

that,  for  this  case,  the  l ine  shape  is  a Lonmtz  curve  (see  for  example, 

Ref.  33).  , 

, . A v'-  ' 

/'■  / v«—  . *•.  ) — v — : .. 

no 


X // 


L (v- -■><)'-  H/SxYaY:-J 


Here.Av  is  the  line  width  (see  sketch) , giver,  by 

Z. 
ir 


A 


r 


(17) 


where  Z is  the  molecular  collision  frequency. 

For  absorption  at  lint?  center,  v-v  - 0 and,  from  eqs.  (Id), 


(1G)  and  (17), 
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Mode  I,  hence  increasing  the  population  inversion.  The  over-all, 
complex  vibrational  relaxation  mechanism  will  be  addressed  in  more  detail 
In  Section  C. 

C.  CALCULATION  OF  GASDYNAMIC  LASER  PERFORMANCE 
1.  Analysis 

In  this  section,  v/e  now  discuss  the  details  of  calculating 

gasdynamic  laser  gain  and  power.  An  accurate  prediction  of  the  population 

inversion  in  an  expanding  gas  encompasses  an  accurate  solution  of  the 

nonequil ibrium,  gas  dynamic  flow  through  a convergent-divergent  nozzle. 

Previous  numerical  solutions  of  high -temperature,  quasi-one-dimensional 

nonequil ibrium  nozzle  flows  where  vibrational  and  chemical  nonequil ibrium 

conditions  prevail  both  upstream  and  downstream  of  the  throat  have  been 

obtained  with  some  effort  by  means  of  a number  of  independent,  steady- 
33 

flow  analyses  . On  the  other  hand,  a useful,  alternative  approach 
has  recently  been  presented  in  References  34  and  35,  namely,  a time- 
dependent  technique  which  entails  the  finite-difference  solution  of  the 
unsteady,  quasi-one-dimensional  conservation  equat’ono  in  steps  of  time, 
starting  from  assumed  initial  distributions  of  the  flow-field  variables 
throughout  the  nozzle.  The  steady-state  nonequilibrium  nozzle  flov/  is 
approached  at  large  times.  The  main  virtue  of  the  time-dependent 
solution  is  its  simplicity;  the  governing  equations  are  solved  by 
means  of  a simple,  explicit,  finite-difference  scheme,  no  additional 
mathematical  methods  are  necessary  to  treat  special  contingencies 
that  can  arise  in  the  analysis  of  nonequilibrium  nozzle  flows,  no 
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difficulties  are  encountered  in  starting  the  nonequilibrium  flow 
from  equilibrium  reservoir  conditions,  very  large  spacings  can  be 
taken  betv/een  grid  points  along  the  nozzle,  and  the  unknown  critical 
mass  flow  is  automatically  approached  at  large  times.  As  a result, 
no  problems  are  encountered  in  the  throat  region  with  the  saddle- 
point  singularity  vcich  plague^  the  steady-flow  formulation. 

The  discussion  that  follows  is  patterned  after  the  technique 
described  in  Reference  8.  Consider  the  grid  point  system  shown  schematically 
in  Figure  9 for  a nozzle  of  specified  shape,  A = A(x),  where  A is  the 
cross-sectional  area.  There  are  two  sets  of  equally  spaced  grid  points; 
one  set  involves  relatively  fine  spacing  from  the  reservoir  to  slightly 
downstream  of  the  throat  in  order  to  enhance  numerical  accuracy  in 
the  flow  region  where  nonequil ibrium  phenomena  develop  at  a fast  rate, 
and  the  second  set  involves  coarse  spacing  further  downstream  where 
nonequilibrium  phenomena  are  usually  progressing  at  a slow  rate.  This 
grid  network  differs  from  that  of  Reference  34,  which  considered  equally 
spaced  grid  points  throughout  the  nozzle.  In  the  present  results  15 
and  5 spaces  were  employed  for  the  fine  and  coarse  grids,  respectively, 
i.e.  a total  of  21  grid  points  were  employed  throughout  the  nozzle. 

The  analysis  begins  by  assuming  initial  values  for  the 
density  p,  velocity  u,  temperature  T,  and  the  nonequilibrium  quantities 
q^  at  all  grid  points  throughout  the  nozzle.  The  precise  quantities 
represented  by  will  be  discussed  later.  Equilibrium  reservoir 
conditions  for  p,  T and  q,.  are  assumed  at  the  first  grid  point  and  field 
fixed,  invariant  with  time.  At  each  internal  grid  point,  spatial 
derivatives  are  then  computed  from  central  finite-differences, 
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^cj.  _ $(x+A/)  -pfx-Ax  ) 

4 x 

s (^X)-2 


(x-Ax) 


(21) 


(22) 


where  g stands  for  p,  u,  T and  q^.  The  nisiibers  for  the  x-derivatives 
obtained  from  equations  (21)  and  (22)  are  then  substituted  into  the 
quasi -one-dimensi onal , unsteady,  conservation  equations 


Continuity 

§f  --  -* 

(23) 

Momentum 

& = - 

(24) 

Energy 

(25) 

Rate 

3^'  = U’A  (f,  Zpj)  - u 

(26) 

State 

(27) 

where  q.  can  be  the  population  or  energy  per  unit  mass  associated 
with  a given  vibrational  energy  level  (or  grouping  of  levels),  and 
w^  is  the  net  internal  rate  °f  production  of  q^  due  to  detailed 
vibration-vibration  and  vibration-translation  collisions  within  the 
gas  mixture. 

Values  of  the  time  derivatives,  |£,  are  directly  obtained 
from  equations  (23)  through  (27),  and  these  numbers  are  subsequently 
inserted  into  an  expression  containing  the  first  three  terms  of  a 
Taylor's  series  expansion  in  time,  i.e., 

CjCJ-s-At)  f ( jj. ^ J.fOjf.  ~z 

Equation  (28)  allows  the  direct  computation  of  all  flow-field 


(28) 
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variables  at  all  Internal  grid  points  at  time  (t+At)  from  the 
known  flow- field  variables  at  time  t.  Values  for  g(t)  are  known; 
as  previously  noted,  values  for  (ff)t  are  obtained  from  equations  (23) 
through  (27).  numbers  for  (^f)t  can  be  obtained  by  differentiation 
of  equations  (23)  through  (26)  with  respect  to  t;  however,  such  an 
operation  introduces  mixed  derivatives.  |j|t.  values  for  which  are 
simply  obtained  by  differentiation  of  equations  (23)  through  (26) 
with  respect  to  x.  Hence,  starting  from  assumed  initial  values  for  the 
flow-field  variables  throughout  the  nozzle,  new  values  of  the  flow-field 
distributions  are  computed  in  steps  of  time  from  equation  (2D), 
continuing  for  many  time  steps  until  the  steady-state  solution 
is  approached  at  large  times  (where  ff  and  *-f  both  approach  zero). 

For  the  present  study,  this  steady-state  solution  is  the  desired 

result. 

The  most  recoct  calculations  using  the  above  approach  employ 
a slightly  modified  technique  discussed  in  References  32  and  35. 

The  modification  consists  of  using 

> + C J+Le  A (29) 

in  lieu  of  eq.  (28).  Here,  the  time  derivative  is  not  evaluated 
at  time  t as  in  eq.  (28);  rather,  an  average  value  between  t and  (t+At) 
is  utilized.  This  average  value  is  obtained  from  the  general  method 
of  HacCormack.  who  has  shown  that  the  general  scheme  is  of  second- 
order  accuracy.  The  results  obtained  using  either  approach  yield  the 
same  results35.  However,  by  using  eq.  (29)  with  MacCormack's 
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finite-difference  scheme,  a considerable  simplification  in  computer 
programing  and  reductions  i.i  execution  time  are  obtained.  See  Kef. 

35  for  more  details. 

The  value  of  at  employed  in  equations  (28)  or  (29)  must 
satisfy  two  stability  criteria*3** 

At  £ ax/  (u+a)  (30) 

where  a is  the  local  frozen  speed  of  sound,  and 

At  < \ (31) 

where  t is  the  characteristic  relaxation  time  for  the  fastest, 
finite-rate  molecular  relaxation  process  occurring  within  the 
mixture.  Equation  (30)  is  the  Courant-Friedrichs-Levy  stability 
criterion,  which  in  physical  terms  states  that  At  must  be  no  greater 
than  the  time  required  for  a sound  wave  to  propagate  from  one  grid  point 
to  the  next,  and  equation  (31)  states  that^t  must  be  no  greater  than 
the  characteristic  time  of  the  finite-rate  molecular  energy  exchange 
processes  occurring  within  the  flowing  gas.  This  latter  stability 
criterion  is  plausible  cri  physical  grounds,  and  has  been  dimonstrated 
empirically  in  Reference  34, 

Because  central  finite-differences  are  employed,  equations 
(21)  through  (28)  can  be  applied  to  the  internal  points  only.  For 
each  time  step,  values  for  p,  u,  T and  q.  at  the  last  grid  point 
(nozzle  exit)  and  u at  the  first  grid  point  (nozzle  inlet)  are  obtained 
by  extrapolation,  as  discussed  in  Reference  34.  Also,  for  the  actual 
numerical  computations,  equations  (21)  through  (30)  are  rewritten 
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In  nondlmcnsional  form  (Reference  34)  so  that  the  magnitudes  of  all 
variables  would  be  approximately  the  same. 

2.  Vibrational  Model 

The  vibrational  model  assumed  in  the  present  study  is  a 
reasonable  approximation  for  the  many  detailed  translation-vibration 
and  vibration-vibration  molecular  collisional  energy  exchanges 
which  can  occur  in  C02-N2-H20  or  C02”^2”He  mixtures.  Figure  10  shows 
a schematic  diagram  of  the  major  participating  vibrational  energy 
levels  for  CO2  and  N2  at  the  temperature  considered  in  the  present 
investigation  (from  300  to  1 500°K) . In  1969,  Taylor  and  Bitterman^8 
made  a thorough  study  and  compilation  of  the  collisional  transition 
probabilities  for  this  system.  Their  results  show  an  extremely  fast, 
resonant,  vibrational  energy  exchange  between  the  v = 1 level  of  lb,  and 
the  (001)  level  of  C02,  as  well  as  a very  fast  exchange  between  the 
(100)  and  (020)  levels  of  C02  due  to  Fermi  resonance.  In  addition, 
vibrational  energy  is  rapidly  transferred  among  the  lower  excited  levels 
of  the  degenerate  mode  v2  in  C02  due  to  the  nearly  equal  spacing  of 
these  levels.  Hence,  these  fast  transitions  appear  to  justify  a 
vibrational  model  which  groups  the  participating  levels  into  two  "modes" 
(Modes  I and  II  in  Figure  10)  which  are  assumed  in  equilibrium  within 
themselves,  but  not  with  each  other. 

More  recent  kinetic  data  obtained  since  1969  are  summarized 
in  References  25  and  26.  In  fact,  Reference  25  is  reproduced  in  its 
entirety  as  Appendices  A and  B in  these  Notes.  The  reader  should 
consult  Appendices  A and  B for  up-to-date  kinetic  rates  before  attempting 


any  calculations  of  his  own. 

In  the  following  discussion  of  the  vibrational  model,  a CUg- 
^-HgO  mixture  will  be  assumed  for  purposes  of  illustration.  COg- 
Ng-He  Is  handled  In  the  same  manner,  with  HgO  replaced  by  He  In  all 
subsequent  equations. 

The  net  vibrational  energies  per  unit  mass  contained  within 

each  "mode",  (ev1't,^  1 and  (evibhl*  are  c^osen  as  dependent  non- 
equilibrium  variables;  hence,  equation  (26)  represents  two  rate 
equations  with  q^  3 ( evib) j and  (eyib)jj  respectively.  The  relaxation 
of  ches'  energies  is  assumed  to  be  described  by  the  simple-harmonic 
oscillator  relaxation  equations; 
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(32a) 


when  q.  3 (evib)It  and 


(32b) 


when  q^  3 ( evib^ 1 1 ' *n  ecluat'’ons  (32a)  and  (32b),  (ey^b)j^  and 
( evib)  1 1 are  the  oclu''^br'ium  vibrational  energies  that  would  be 
contained  in  Modes  I and  II  respectively  at  the  local  gas  trans- 
lational temperature  T,  and  Tj  and  t j j are  the  characteristic  relaxation 
times  for  Modes  I and  II  respectively,  tj  and  t j j are  averages  which 
characterize  the  net  rate  of  energy  transfer  into  and  out  of  Modes  I 
and  II;  this  energy  transfer  is  assumed  to  be  governed  by  the  major 
transitions  (heavy  arrows  in  Figure  10)  which  are  identified  with  the 
relaxation  times  t , t.  , t_.  These  relaxation  times  are  themselves  averages 
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of  the  detailed  C02-C02,  C02-N2,  C02-H20,  N2«N2  and  N2«H20  collisions; 
such  averages  for  a mixture  of  gases  are  obtained  from  the  "parallel 
resistance"  mixture  rule: 
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Xj  denotes  the  mole  fraction  of  chemical  species  i.  For  the  present 
results,  values  for  the  denominators  of  the  terns  on  the  r.h.s.  of 
equations  (33a,  b and  c)  are  obtained  directly  from  Reference  38. 
(However,  the  reader  should  consult  Appendices  A and  B of  the  present 
notes  for  more  recent  kinetic  rates.)  In  turn,  the  average  relaxation 
times  for  Modes  I and  II  are  obtained  from 


(34) 


+ 


(35) 


These  values  are  subsequently  employed  in  equations  (32). 

Parenthetically,  the  general  quantity  t can  be  interpreted 
as  a mean  time  required  for  a single  particle  to  make  a transition 
from  one  state  to  another  state  due  to  collisions  with  other  particles. 
In  turn,  can  be  interpreted  as  the  number  of  transitions  per  second 
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per  particle.  Hence,  on  a physical  basis,  equations  (33a,  b and  c) 
simply  state  that  the  total  number  of  transitions  per  second  per  particle 
due  to  collisions  with  all  species  present  in  the  mixture  is  the  sum 
of  the  transitions  per  second  per  particle  due  to  collisions  with  each 
individual  specie-  Equations  (34)  and  (35)  have  similar  physical 
interpretations,  except  now  the  "transitions"  are  from  one  mode  to 
another. 

The  above  model  for  treating  the  complex  vibrational 
energy  transfer  processes  Is  approximate,  and  its  limitations  are 
discussed  in  Reference  11. 

It  is  sufficient  to  state  here  that  the  model  is  intended 
only  for  the  calculation  of  population  inversions  in  COg-i^-^O  or  He 
mixtures;  it  is  not  necessarily  valid  for  other  gases,  nor  can  it 
be  used  when  substantial  amounts  of  radiative  energy  interact  with 
the  gas. 

An  alternative  to  the  above  model  is  the  exact  solution 

39 

of  the  master  relaxation  equations  for  each  energy  level  . However, 
for  many  cases  of  interest,  comparisons  with  experiment  (to  be  discussed) 
show  that  such  a detailed  treatment  is  not  absolutely  necessary,  and 
that  the  present  model  yields  reasonable  accuracy. 

Referring  again  to  equations  (32a)  and  (32b),  (ey^)j  and 
( evib ) 1 1 are  treated  as  the  dependent  nonequilibrium  variables.  From 
these  energies,  different  vibrational  temperatures,  (Ty^b)j  and 
(TVib) i i * are  defined  from  the  following  equilibrium  relations: 
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In  equations  (36a  and  b),  cCQ  and  are  the  mass  fraction  and 
specific  gas  constant  repsectlve'y  for  C02 ; cN  and  are  defined 
similarly  for  N«;  v2  and  v3  are  the  normal  vibrational  frequencies 
of  the  symmetric  stretching,  the  bending,  and  the  asymmetric  stretching 
modes  respectively  for  C02;  andv'is  the  normal  vibrational  frequency 
for  N2.  In  turn,  (Tvib)I  and  (Tvib)n  are  used  to  compute  populations 
of  energy  levels  within  Modes  I and  II  assuming  a Boltzmann  distribution 
locally  within  each  mode.  For  example,  the  population  of  the  (001) 
level  in  C02  is  obtained  from 
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where 


<?=£/-  i^,u- Tyrt*]  £/-£~a^u h*J 


(30) 


and  the  population  of  the  (100)  level  in  C02  is  obtained  from 
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In  the  above  formulation,  the  role  of  H20  or  He  as  a catylist 
enters  the  analysis  through  its  presence  In  equations  (33a,  b and  c). 

As  a final  note  with  regard  to  the  present  analysis,  (eyib)j 
and  (evi-b)jj  are  inherent  parts  of  the  time-dependent  nonequilibrium 
nozzle  flow  solution,  and  are  computed  at  each  time  step  during  the 
transient  approach  to  steady-state  conditions.  Then,  after  the 
steady-state  is  achieved,  the  vibrational  temperatures  and  populations 
are  computed  from  equations  (36)  through  (39). 

3.  Results 

The  numerical  nature  of  the  time-dependent,  vibrational 
nonequl ibrium  nozzle  solution  is  clearly  shown  in  Figure  11,  which 
Illustrates  the  shape  of  the  (eyjb)jj  distribution  through  a wedge 
nozzle  at  several  times  during  its  approach  to  the  steady  state.  In 
Figure  11,  R is  the  specific  gas  constant  for  the  mixture.  The 
dotted  line  shows  the  assumed  initial  distribution  at  t = 0,  the 
solid  lines  show  distributions  after  400,  800  and  1000  time  steps, 
and  the  solid  dots  show  the  final,  steady-state,  nonequilibrium 
distribution  obtained  after  approximately  3000  time  steps.  Similar 
transient  variations  and  the  subsequent  approach  to  the  steady-state 
occur  for  all  other  flow  variables  throughout  the  nozzle. 

The  final,  steady-state  (eyib)jj  distribution  shown  in 
Figure  11,  taken  in  conjunction  with  the  simultaneous  steady-state 
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(e  1b),  distribution,  loads  to  the  steady-flow  result,  shown  in 
Figures  12  through  14.  In  Figure  12.  the  two  vibrational  temperature 
distributions  are  compared  with  the  translational  temperature  distribution. 
Note  that  Mode  11  relaxes  much  slower  than  Mode  I.  and  that  both 
Modes  I and  11  are  not  in  equilibrium  with  the  local  translational 
temperature  except  near  the  nozzle  reservoir.  Also,  note  that 
(T  begins  to  diverge  from  T in  the  subsonic  section  of  the 

nozzle,  with  [(Tv1b)n-T]  - 1*>’K  at  the  throat.  Hence,  nonequilibrium 
effects  upstream  of  the  throat  are  noticeable.  In  Figure  13,  the 
populations  N001  and  N,00  of  the  (001)  and  (100)  levels  respectively 
in  C02  are  shown  as  functions  of  distance  along  the  wedge  nozzle. 

At  approximately  0.26  cm  downstream  of  the  throat,  the  populations 
Of  the  (001)  and  (100)  levels  are  equal,  and  a substantial  populate 
inversion  develops  further  downstream.  This  population  inversion 

is  Shown  in  Figure  14  in  terms  of  (!,00i -;,ioo) /l!C02  “nere  ,s  th’ 
incal  number  density  of  C02  particles.  If  the  nozzle  exhausts  into 

a constant  area  duct  as  shown  in  Figure  14.  tne  inversion  will 
continue  to  increase  beyond  the  nozzle  exit  until  Mode  1 equilibrates 
with  T.  beyond  which  the  inversion  will  slowly  decrease  due  to  the 
moderate  deactivation  of  Mode  II.  As  emphasized  in  References  7 
and  38.  the  equilibration  of  Mode  I with  T is  enhanced  by  the  presence 

of  H20  in  the  mixture. 

4.  Conparison  with  Basov 

Basov  et  al6  have  published  results  for  population  inversions 
created  in  . supersonic  expansion  of  a C02-N2  mixture  through  a hyporboli 


THROUGHOUT  NOZZLE 


Fig.  14  STEADY-STATE  DISTRIBUTION  OF  POPULATION 
INVERSION  THROUGHOUT  NOZZLE 


nozzle;  In  contrast  to  the  present  results  shown  in  Figures  11 
through  14,  Basov  did  not  include  the  effect  of  HgO  on  the  population 
inversions.  Also  in  contrast  to  the  present  analysis,  Basov  assumed 
equilibrium  conditions  from  the  reservoir  to  the  throat.  In  order  to 
provide  a comparison  with  the  results  of  Reference®,  results  are 
obtained  with  the  present  time-dependent  analysis  for  a nozzle  shape, 
p , T0  and  mixture  ratio  identical  to  those  of  Reference®.  Figures  15 
and  16  show  the  resulting  comparison;  Figure  15  illustrates  the 
populations  N^qq  and  NQg1  as  a function  of  distance  through  the  nozzle, 
and  Figure  16  shows  the  population  inversion  distribution  throughout 
the  nozzle.  NQ  is  the  total  reservoir  number  density  for  the  gas 
mixture.  A study  of  Figure  15  shows  that  Basov  predicts  a slower 
relaxation  of  the  (001)  level  and  a faster  relaxation  of  the  (100) 
level  in  comparison  to  the  present  results.  Consequently,  as  seen 
in  Figure  16,  Basov's  results  show  a population  inversion,  whereas  for 
the  given  conditions  the  present  results  predict  no  population  inversion. 
The  present  results  embody  detailed  information  on  transition 
probabilities  reported  in  Reference  38,  whereas  Basov  used  some- 
what cruder  information  for  transition  probabilities.  A comparison 
of  these  two  sets  of  transition  data  show  approximate  agreement  for 
the  collision  probability  for  deactivation  of  the  (001)  level;  on 
che  other  hand,  Basov  used  prcbabilities  for  deactivation  of  the 
(010)  level  (hence  the  (100)  level,  which  is  closely  coupled  by  fast 
resonant  exchanges)  which  are  more  than  an  order  of  magnitude  larger 
than  those  presented  in  Reference  38.  Consequently,  the  differences 
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shown  In  Figures  15  and  16  can  be  attributed  mainly  to  differences 
In  relaxation  times  used  for  the  two  analyses,  Taylor  and  Bltterman's 
thorough  compilation  of  transition  probabilities  (Refs.  7 and  38)  is 
believed  to  contain  the  most  reliable  data  currently  available  in 
the  literature.  These  data  have  been  very  recently  up-dated  by 
Anderson  In  References  25  and  26.  (See  Appendices  A and  ITof  the 
present  Notes).  Also,  these  results  clearly  demonstrate  that 
population  inversions  in  COg-Ng  mixtures  without  a catalyst  are 
difficult  to  obtain.  Indeed,  In  subsequent  sections  of  these  notes, 
we  will  see  that  the  presence  of  H20  or  He  is  vitally  necessary  for 
the  attainment  of  reasonable  population  inversions. 

As  a point  of  interest,  Basov's  assumption  (Ref.  6)  of 
equilibrium  conditions  at  the  nozzle  throat  is  reasonable  for  the 
nozzle  shape  he  considered;  for  this  nozzle,  the  subsonic  approach 
to  the  throat  is  relatively  long  and  smooth,  thus  enhancing 
equilibrium  conditions.  In  fact,  the  present  results,  v/hich  include 
nonequilibrium  conditions  upstream  of  the  throat,  indicate  that  C(Ty^b) j i -T] 
is  28°K  at  the  throat  for  Basov's  conditions.  This  should  be  compared 
with  the  120°K  difference  obtained  for  the  shorter  nozzle  considered 
in  Figure  11,  where  nonequilibrium  effects  upstream  of  the  throat 
are  considerable. 

5.  Coupled  and  Uncoupled  Flows 

With  regard  to  calculations  of  population  inversions  in  an 
expanding  nozzle  flow,  there  has  been  some  question  regarding  the 
error  induced  by  uncoupling  the  gas  dynamic  flov/  field  from  the 
ibrational  relaxation  phenomena.  For  example,  the  calculation  of 
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population  inversions  would  be  greatly  simplified  if  the  pressure, 
temperature  and  velocity  distributions  were  obtained  from  a 
calculation  of  frozen  or  equilibrium  nozzle  flows,  and  then  used 
to  solve  independently  the  vibrational  rate  equations.  The  flow 
variable  most  severely  affected  by  these  uncoupled  solutions  is 
translational  (static)  temperature.  Since  the  relaxation  rates 
depend  exponentially  on  T,  and  because  tries  to  equilibrate 
with  T in  the  supersonic  stream,  the  theoretical  population  inversions 
obtained  from  such  uncoupled  solutions  may  contain  an  inordinate  error. 

In  order  to  examine  this  question,  a numerical  experiment  is 
performed  using  the  present  time-dependent  analysis  to  compute 
population  inversions  for  the  following  cases:  (1)  Coupled,  where 

the  finite-rate,  vibrational  energy  relaxation  is  fully  included  in 
the  gas  dynamic  energy  equation,  (2)  uncoupled  (equilibrium),  where 
p,  T and  u are  obtained  from  an  expansion  assuming  infinite-rate, 
vibrational  equilibrium,  and  (3)  uncoupled  (frozen),  where  p,  T and 
u are  obtained  from  an  expansion  assuming  a zero  vibrational  rate,  i.e., 
frozen  flow.  In  case  (3),  two  sets  of  results  are  obtained,  where  the 
ratio  of  specific  heats  y equals  1.4  and  1.3.  The  value  of  1.4  properly 
corresponds  to  a vibrationally  frozen  expansion  of  CO2  and  whereas 
the  value  of  1.3  is  of  interest  only  because  it  is  equal  to  the 
equilibrium  value  of  the  ratio  of  specific  heats  for  the  gas  mixture 
in  the  reservoir.  For  all  of  the  above  cases,  a simple  wedge  nozzle 
was  assumed,  followed  by  a constant  area  duct.  This  geometry  and  the 
reservoir  conditions  are  described  in  Figure  17.  Note  that  this  nozzle 


has  the  same  area  ratio  but  twice  the  length  of  the  wedge  nozzle  previously 
considered  in  Figures  11  through  14. 

Results  from  this  numerical  experiment  are  shown  in  Figures  17 
through  20.  Figure  17  contrasts  the  translational  temperature  distri- 
butions obtained  for  the  different  assumptions  of  infinite,  finite  and 
frozen  rates;  such  a comparison  is  classical,  and  is  presented  in 
order  to  help  interpret  the  subsequent  results.  Also,  Figure  17  shows 
that,  for  the  present  conditions,  the  actual  nonequilibrium  expansion 
is  closer  to  equilibrium  than  frozen.  In  turn,  a comparison  of  the 
curves  shown  in  Figure  18  indicates  that  (Ty^b)j  is  more  strongly  affected 
than  (Tv^)jp  Indeed,  the  upper  level  appears  to  be  insensitive  to 
the  differences  in  translational  temperature  between  the  different  cases 
despite  the  exponential  variation  of  relaxation  time  with  T.  The  reason 
for  the  differences  in  (Tvib)j  is  the  strong  tendency  of  the  lower  level 
to  equilibrate  with  T,  which  is  different  for  each  of  the  cases. 

The  differences  in  (Ty^)j  shown  in  Figure  18  are  magnified 
when  [J.|qq  is  computed,  as  shown  in  Figure  19.  Also,  as  expected 
from  the  very  small  differences  in  only  small  differences 

occur  in  . 

The  differences  shown  in  Figure  19  translate  into  the 
differences  in  population  inversion  shown  in  Figure  20.  Among  the 
simplified,  uncoupled  cases.  Figure  20  illustrates  that  the  p,T  and  u 
variations  obtained  from  an  equilibrium  flow  calculation  provide  the 
closest  comparison  with  the  coupled  case.  However,  there  is  still  a 
noticeable  error  induced  by  all  the  uncoupled  cases.  With  regard  to 
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the  frozen  cases,  one  might  be  able  to  simulate  the  coupled  results  by 
"proper"  choice  of  an  effective  y between  1.3  and  1.4,  however,  such 
a value  of  y is  not  known  a priori,  and  its  proper  choice  would  be 
strictly  fortuitous  without  first  knowing  the  anser. 

Reflecting  upon  the  results  shown  in  Figures  17  through  20, 

i 

the  accurate  calculation  of  population  Inversions  in  an  expanding 
gas  apparently  requires  an  analysis  which  fully  takes  into  account 
the  coupling  between  the  finite-rate  nonequilibrium  processes  and 
the  gas  dynamic  flow. 

It  is  interesting  to  note  that,  for  the  same  reservoir 
conditions,  the  peak  population  Inversion  shown  in  Figure  14  is  23 
percent  higher  than  the  peak  shown  in  Figure  20  for  the  coupled 
case.  The  larger  inversion  shov/n  in  Figure  14  is  due  solely  to 
the  shorter  nozzle  length,  and  is  to  be  expected  due  to  a greater  degree 
of  nonequilibrium  in  the  shorter  nozzle  flow. 

6.  Recapitulation 

The  previous  sections  have  described  the  time-dependent  technique 
for  calculating  population  inversions  in  an  expanding  gas,  and  have 
given  some  typical  results  obtained  with  an  approximate  vibrational  model. 
The  kinetic  rates  employed  in  generating  the  results  shown  in  Figures  11 
through  20  are  given  in  Appendix  A of  Reference  40,  which  contains  simple 
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correlations  of  the  rates  from  Taylor  and  Bitterman  . However,  these 
correlations  have  been  improved  as  discussed  in  References  25  and  26, 
and  the  reader  is  again  urged  to  consult  Appendices  A and  B of  the 
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present  Notes  for  up-dated  kinetic  rates. 

Another  detailed  numerical  solution  for  GDI's  can  be  found 
In  Reference  41,  where  Munjee  has  Incorporated  a more  detailed 
kinetic  model  than  used  In  Reference  8.  Munjee 's  kinetic  model  is 
patterned  after  that  of  Basov  et  al^-,  however,  unlike  Basov,  Munjee  in- 
cludes the  effect  of  a catalyst  (HgO  or  He)  in  the  master  rate 

equations.  Munjee's  work  serves  as  the  basis  for  the  theoretical 
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results  quoted  by  Christiansen,  Tsongas  and  Buonadonna  ’ . In 

addition  to  these  detailed  numerical  solutions,  a useful  approximate 
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analysis  has  been  developed  by  Hoffman  and  Vlases  . Also,  recent 
results  using  the  sudden-freezing  approximation  have  been  reported  in 
Reference  28. 

7,  Power  Extraction 

In  a gasdynamic  laser,  the  energy  avilable  for  laser 
power  extraction  is  contained  in  the  vibrational  energies  of  the  excited 
N2  and  C02(v3)  at  the  exit  of  the  supersonic  nozzle;  the  challenge 
is  to  optically  extract  as  much  of  this  available  energy  as  possible 
before  it  deactivates  and  "leaks"  into  translational  and  rotational 
energy  due  to  molecular  collisions.  Moreover,  the  extraction  process 
should  also  produce  a good  quality,  near-diffraction-limited  beam. 
Therefore,  the  optical  design  of  the  laser  cavity  is  of  utmost  importance. 
Some  considerations  dealing  with  the  optical  cavity  will  be  discussed 
in  Section  III  of  these  Notes,  and  some  of  the  aerodynamic  distrubances 
which  Influence  cavity  design  and  beam  quality  will  also  be  treated  in 
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Section  111.  However,  let  us  now  consider  power  extraction  from  a 
purely  energy  standpoint. 

The  maximum  available  laser  energy  can  be  calculated  as 
follows.  Consider  a point  In  the  nozzle  of  a gasdynamic  laser. 

Using  the  notation  of  Reference  8,  the  translational  and  rotational 
temperature  at  this  point  is  T,  the  sum  of  the  vibrational  energy 
contained  in  the  excited  N ^ and  COg  (v^)  is  ey^b  with  an  attendant 
vibrational  temperature  T..,h  , the  sum  of  the  vibrational  energy  in 

VIDji 

the  COgfvj)  and  COgfvg)  modes  Is  ey-b  with  an  attendant  Tyib  , and 
the  population  densities  of  the  upper  and  lower  laser  levels  are  Nqqj 
and  N100  respectively.  In  GDL  flows  in  general,  Tyjb  >Tvibf”T* 

When  an  inversion  exists,  then  by  definition  ilQ01  - N1Q0>0.  The 
population  densities  are  given  by 

/-A 


Mpo/  - /cp')  C 


--  (A/c^/<p)  e~C~'“MTv;b* 


where  Q is  the  partition  function.  We  ask  the  following  question:  If 

energy  is  drained  from  ey-b  , and  ey^b  is  held  constant,  at  v/hat 
value  of  Ty^b  will  the  population  inversion  go  to  zero?  Denote 
this  value  of  T by  (T  ih°  ).  Then, 

VIDji  vidu 


Ho/  ~ Heo  - O - (A/Ci ,A/tp)  Cc 
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Hence, 


j 

) 

) 

1 


^CO/  / - ^/£>c>  /•&  rvtbx 

C ’Tr/'Stf  ) ” 'T''t'hjr  ( ^co/  ) 

When  T > (t  ..°  ),  an  inversion  exists  and  laser  power  can  in 
wnen  vibjj  ' vibjj 

principle  be  extracted.  When  TVibn  < ^vibj^’  n0  inversion  is 
present,  and  no  power  can  be  extracted. 

In  order  to  consider  the  maximum  laser  energy  available,  we 
assume  that  Tv(b  -T.  Then.  (T^M  - T(e001/e,n0>.  Using  this  value 
of  (T  . • ).  we  define  a maximum  available  laser  energy  as 

w ' vibjj  ’ 

~ L T,rtba  ) ~ C Xi+*)1 

where  the  factor  0.409  is  the  quantum  efficiency  for  the  C02  laser  transition 
at  1 0.6^. 

The  quantity  e , is  a convenient  index  to  gage  the  amount 

of  power  that  might  be  extracted  from  a gasdynamic  laser.  Howevar,  in 

reality  the  actual  power  extraction  is  usually  less  than  e^  due  to 

losses  in  the  laser  cavity.  Also,  all  values  of  emax  given  in  the  present 

Notes  are  local  values  at  the  nozzle  exit,  where  is  obtained 

from  the  coupled  analysis  of  Reference  8.  Therefore,  the  present  values 

of  e account  for  the  kinetic  deactivation  losses  in  the  nozzle, 
max 

This  is  in  contrast  to  previous  simpler  but  less  realistic  definitions 
of  maximum  available  power1,  which  have  been  based  strictly  on  the 
vibrational  energy  in  the  reservoir. 


4G 


m 


Emphasis  Is  again  made  that  ema„  is  simply  a convenient 
index  that  represents  an  upper  limit,  and  that  the  actual  power 
extracted  may  be  quite  different.  An  accurate  calculation  of  power 
extraction  must  combine  a detailed  numerical  solution  of  the  non- 
equilibrium supersonic  flow  in  the  cavity  coupled  with  a point-by- 
point  numerical  integration  of  the  radiative  transfer  equation 
C/Tr/Tr  = <? 

in  a direction  normal  to  the  flow.  Here,  6 is  the  actual  gain 

{ proportional  to  the  actual  population  Inversion)  which  takes  into 

account  the  local  depopulation  of  the  H-  and  COg  (v^)  levels,  and  the 

population  of  the  COg  ( ) levels  due  to  interaction  with  the  laser  radiation. 

6 is  less  than  GQ,  which  is  the  small-signal  gain  coefficient  discussed 

previously.  This  calculation  should  be  at  least  a "two-dimensional" 

analysis,  where  I y- varies  In  the  flow  direction  as  well  as  in  the 

beam  direction  perpendicular  to  the  flow.  Such  detailed  calculations 
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have  been  performed  by  Heiche  and  Harris  at  the  Naval  Ordnance 

Laboratory  (NOL),  but  they  are  not  generally  available  in  the  literature. 

However,  for  the  sake  of  comparison,  the  following  example  is  given  to 

compare  e with  the  calculated  actual  power  extracted  from  a GDL  cavity, 
max 

More  details  on  this  example  can  be  found  In  Reference  22. 

Consider  a GDL  cavity  which  is  50.5  cm  wide  in  the  beam  direction 
(transverse  to  the  flow)  and  1.5  cm  high.  This  size  is  representative 
of  the  NOL  3-Hegawatt  Arc  Tunnel  after  it  was  converted  for  GDL  experiments, 
as  described  in  Refereces  11  and  12.  Consider  also  a very  hypothetical 
Master  Oscillator  Power  Amplifier  (HOPA)  arrangement  (see  section  III  of 
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these  Notes)  consisting  of  a 1.5  x 10  cm  rectangular  beam  traversing 
the  cavity  with  five  side-by-side  passes  beginning  at  the  nozzle  exit. 
Hence,  power  is  extracted  over  a cavity  length  of  50  cm.  Assume  that 
the  input  laser  beam  has  a power  P1n  - 10  KW  with  a Guassian  intensity 
distribution  in  the  flow  direction.  Let  PQut  be  the  output  power  after 
the  last  pass.  The  power  extracted  from  the  gas  is  Pfi  » PQUt  - p1n* 
With  this  arrangement,  the  coupled  power  extraction  - nonequilibrium 
gasdynamlc  analysis  of  Helche  and  Harris43  yields  the  results  shown  In 

Table  I for  Pg. 


TABLE  I 


XH20 

emax 

(KJ/lbm) 

at 

nozzle 

exit 

Pe/ni 

(KJ/lbm) 

Percent 

ma?f 

available 

0.01 

23.4 

12.2 

521 

0.035 

19.8 

10.5 

531 

0.07 

15.7 

6.67 

43« 

jh  * 1.5  lb/sec 


For  these  results,  the  reservo  - temperature  and  pressure  were  1800‘K 
and  37.5  atm  respectively,  the  mole  fraction  of  C02  was  0.07,  the  nozzle 
area  ratio  was  Ag/A*  ■ 50,  and  the  nozzle  throat  height  was  h*-0.3  mm. 
Three  cases  with  different  H20  mole  fractions  are  shown  In  Table  I. 

These  results  show  that,  for  a GOl  the  size  of  the  HOL  3-llegawatt 
Arc  Tunnel,  Pe  Is  on  the  order  of  1/2  the  maximum  available  power 
at  the  nozzle'exlt.  This  serves  to  Illustrate  the  difference  between 

43 


iid.. 


HM 


IMWUMJWII 


■W*wp,, , I > yapfs  IBT^W-  ■» 


emax  and  the  actua^  P^er  extracted. 

As  a final  note  on  power  extraction,  several  approximate 
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methods  * are  available  for  predicting  power  output,  generally 
based  on  the  estimation  of  saturation  intensity  within  the  cavity 
(see  section  III).  However,  for  accurate  results  necessary  for  reasonable 
comparison  between  theory  and  experiment,  and  for  detailed  GOL  design, 

t 

a coupled  numerical  solution  such  as  described  above  is  required. 

D.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 


The  previous  discussions  have  dealt  with  calculations  of 
GDL  performance.  Such  calculations  clearly  show  that  population 
inversions,  laser  gain,  and  laser  pov/er  extraction  can  occur  in 
supersonic  expansions  of  C02-N2-H20  or  He  mixtures.  These  calculations 
have  been  generally  confirmed  by  experiment,  and  such  experiments 
are  the  suLj.-t.  of  this  section.  Also,  the  results  shown  in  this 
section  draw  heavily  on  References  11  and  45. 

1.  Kinetic  Rate  Data 

The  vibrational  relaxation  times  used  to  ootain  the  theory 
curves  in  this  section  are  correlations  of  the  data  compiled  in 

reference  38,  with  the  single  exception  of  the  water  rate,  (tc)Cq  q. 
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Based  on  recent  calculation  by  Sharma  the  present  model  assumes  that 

( tcp)co  -H  0 constant  f°r  T ~ 600°K.  This  approximation  is  not 

good  for  T<200°K,  but  such  low  temperatures  are  usually  not  encountered 
in  the  flows  of  interest  here. 

2.  Honeguilibrium  Gas  Dynamics 

In  the  analysis  of  reference  8,  the  kinetic  rate  equations 
are  fully  coupled  with  the  governing  flow  equations  of  continuity, 
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momentum  and  energy.  Nonequilibrium  conditions  are  included  both 
upstream  and  downstream  of  the  nozzle  throat.  The  governing  equations, 
along  with  their  time-dependent  numerical  solution  have  been  described 
In  Section  C of  these  Notes.  Hence,  no  further  elaboration  will  be 
given  here. 

3.  Small -Signal  Gain 

The  fundamental  experimental  measurement  discussed  in  this 
Section  is  the  small -signal  gain  coefficient,  GQ,  defined  as  diyi^* 

GQdx.  Here,  Iv  is  the  intensity  of  a low  power  C02  diagnostic 
laser  beam  at  10.6  pm  incident  on  an  element  of  gas  of  thickness  (dx) 
in  the  supersonic  flow,  as  shown  schematically  in  Fig.  21.  In  principle, 
the  intensity  increase,  dlv,  yields  a measurement  of  GQ.  The  intensity, 
Iv,  must  be  small  enough  such  that  the  radiation  field  does  not  disturb 
the  population  inversion,  i.e.,  such  that  (!^ooi"'hoo^  co^^sion 
dominated.  In  this  case,  GQ  is  directly  proportional  to  the  population 
inversion.  A detailed  discussion  and  derivation  of  the  small -signal 
gain  equation  is  given  in  Section  C.  With  the  appropriate  rotational 
constants,  a close  approximation  of  eq.  (20)  for  the  P(20)  transition 
of  the  10.6  band  is 


r - 

* ^ 3,  vc 


^ “ tyoo  ) ( T ) ^ (40) 


•6 

where  A “ 10.6  x 10  m,  *s  the  spontaneous  radiative  lifetime  * 5.38 
sec,  v is  the  collision  frequency,  and  T is  in  °K.  Equation  (40) 

V# 

assumes  pressure  line  broadening  only;  Doppler  broadening  is  negligible 
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for  the  pressures  of  Interest  at  the  nozzle  exit  ( 50  torr).  The 
quantities  In  Eq.  (40)  are  not  limited  to  a specific  set  of  units; 


1 however,  they  must  be  consistent,  such  as  the  SI  system.  The  numerical 

fU 

f constants  in  the  last  two  factors  have  units  of  °K. 


Thus,  theoretical  vibrational  population  Inversions  (Njor^lOO^ 
are  obtained  from  the  analysis  of  reference  8;  in  turn,  these  are  used 
In  Eq.  (40)  to  give  theoretical  predictions  of  Gq  for  the  P( 20)  transitions 
of  the  10. 6u  band.  For  other  P(J)  transitions  In  the  same  band, 

Eq.  (40)  must  be  appropriately  ratloed,  assuming  rotational  equilibrium. 

The  purpose  of  the  remainder  of  this  section  Is  to  compare  the  theoretically 
predicted  GQ  with  experimental  measurements. 

4 . NOL  3-Megawatt  Arc  Tunnel 

The  NOL  3-Megawatt  Arc  Tunnel  Is  an  arc  heated  high  temperature 

12 

supersonic  wind  tunnel  facility  . The  artist's  conception  of  Figure  22 
shows  the  major  conventional  components  of  the  wind  tunnel.  Electrical 
power  Is  delivered  from  an  AC  power  supply  to  the  4 ring,  3 phase  AC  arc 
heater.  Nitrogen  (or  air)  Is  heated  to  the  desired  pressure  and  temperature 
by  the  arc.  In  Its  conventional  wind  tunnel  configuration,  the  high 
temperature  gas  passes  through  the  supersonic  nozzle  into  the  test  cell, 
where  aerodynamic  and  heat  transfer  data  are  usually  obtained.  The  gas 
then  passes  through  the  diffuser,  an  after-cooler,  and  to  a vacuum  pumping 
pi  ant. 
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The  arc  heater  is  of  the  4-ring  3 phase  AC  type  developed 
at  H0L  (See  Figure  23).  Water  cooled  copper  electrodes  and  liners 
are  placed  within  the  steel  pressure  shell.  The  gas  to  be  heated 
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AFTER-COOLER 


Fig.  22  THREE  MEGAWATT  ARC  TUNNEL 


Is  Injected  Into  the  chamber  through  a series  of  holes  around 
one  leg  of  each  electrode.  The  gas  Is  heated  by  the  two  arc 
coluwis  which  (due  to  their  self-induced  magnetic  field)  rotate  around 
the  electrodes.  The  arc  heater  is  capable  of  operation  over  a range 
of  temperatures  from  1500°K  to  6000°K  and  pressure  of  10  to  70 
atmospheres.  Extensive  operating  experience  at  pressures  of  10  to  35 
atmospheres  give  heater  efficiencies  of  25  to  60  percent  depending  upon 
the  temperature  and  mass  flow  required. 

For  the  small-signal  gain  measurements,  the  conventional 
arrangement  shown  In  Figure  22  is  modified.  Specifically,  a mixing 
chamber,  supersonic  rapid  expansion  nozzle  and  constant  area  duct 
are  substituted  for  the  conventional  wind  tunnel  nozzle.  This  modified 
arrangement  is  shown  In  the  schematic  of  Figure  24  and  the  photo  of 
Figure  25.  In  the  mixing  chamber,  CO2  and  liquid  1^0  are  injected  and 
mixed  with  the  hot  from  the  arc  heater.  The  hot  gas  mixture  passes 
through  the  supersonic  rapid  expansion  nozzle,  where  the  population 
Inversion  is  created.  The  two-dimensional,  contoured,  minimum  length 
nozzle  has  a throat  height  of  1 mm,  an  inviscid-core  area  ratio  of  20, 
and  a length  transverse  to  the  flow  of  585  mm.  The  flow  passes  through 
an  essentially  constant  area  duct  which  is  provided  with  a series  of 
7 viewing  ports  along  its  length.  The  duct  is  slightly  divergent  to 
account  for  viscous  effects.  The  output  beam  (10.6  pm,  P(20)  transition) 
from  a homemade  diannostic  COg  gas  laser  (see  Figure  26)  is  directed 
into  the  duct  through  an  IRTRAN  window,  reflected  from  a gold-coated 
mirror  on  the  opposite  side  of  the  duct,  and  passed  back  out  the  window 
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Fin.  24  TOP  VIEW  SCHEMATIC  OP  THE  NOL  3-MEGAWATT  ARC  TUNNEL  MODIFIED  FOR  THE 
PRESENT  EXPERIMENTS 
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Into  a broadband  cw  laser  power  meter  (Coherent  Radiation  Labs 
Model  201).  The  small  signal  gain  along  the  duct  may  be  measured  by 
moving  the  window  and  mirrors  so  that  the  diagnostic  laser  beam 
passes  through  the  desired  point  location. 

The  reservoir  temperature  of  the  gas  mixture  entering  the 

supersonic  nozzle  Is  determined  by  use  of  a one-dimensional  mass 

flow  relationship  from  a sonic  throat  condition.  The  measured  arc 

heater  pressure  and  mass  flows  of  the  individual  gases  (N2,  C02,  and 

H20)  are  used  with  specific  heat  ratios,  compressibility  factors,  and 

gas  constants  for  the  mixture  to  calculate  the  nozzle  reservoir 

temperature.  The  procedure  is  similar  to  that  normally  used  in  high 

48 

temperature  wind  tunnel  systems  . 

The  spectral  output  of  the  diagnostic  ".02  gas  laser  was 
measured  with  the  aid  of  a Jarrell -Ash  model  (8200D)  spectrometer 
and  a liquid  nitrogen  cooled  detector.  Trie  laser  was  found  to  operate 
predominately  on  the  P(2Q)  transition;  the  P( 18)  and  P(lu)  transitions 
also  occurred,  but  were  weak  enough  to  be  ignored  in  the  small -signal 
gain  measurements  and  in  the  analytical  calculations. 

The  sequence  of  events  leading  to  a small -signal  gain 
measurement  is  as  follows.  First,  the  diagnostic  laser  and  optical 
system  are  aligned  to  give  the  maximum  power  Incident  on  the  power  meter. 
Then,  the  arc  heater  is  turned  on,  heating  pure  fl2  at  a temperature 
somewhat  above  the  final  reservoir  temperature.  When  steady  flow 
conditions  are  reached  in  the  heater,  cold  C02  and  liquid  li20  are 
Injected  into  the  mixing  chamber  (see  Figures  24  and  25)#  The  power 
meter  continuously  records  the  power  level  of  the  beam  after  it  has 
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doubly  traversed  the  duct.  The  beam  power  measured  during  the 
steady  flow  of  the  N2-C02-H20  mixture  and  the  power  WQ  measured 
during  the  preliminary  flow  of  pure  N2  yield  GQ  from  the  relationship 
tyWo  * exp  (GqL)  where  L is  the  geometric  path  length  inside  the  duct. 

5.  NOL  12.7  cm  Shock  Tunnel 

The  arc  tunnel  discussed  above  is  one  means  of  providing  a 
reservoir  of  hot,  vibrationally  excited  gas  for  subsequent  expansion 
through  a nozzle.  Another  means  is  to  utilize  the  shock-heated  gas 
behind  a reflected  shock  wave  in  a shock  tube.  In  the  present  experi- 
ments, a conventional  pressure  driven  shock  tube  is  used  to  provide 
reservoir  conditions  behind  a reflected  shock  wave  in  a mixture  of 
C02-N2-He.  A nonequilibrium  population  inversion  is  created  when  this 
shock-heated  mixture  expands  through  a nozzle  mounted  at  the  end  of 
the  tube.  The  flow  then  passes  through  a slightly  diverging  cavity 
where  measurements  of  GQ  are  made,  and  finally  exhausts  into  an 
evacuated  dump  tank.  Some  features  of  the  experimental  apparatus  are 
given  below,  and  a schematic  of  the  shock  tube  arrangement  is  given 
in  Figure  28. 

The  shock  tube  is  12.7  cm  in  diameter  and  utilizes  a single 
diaphragm  to  separate  the  driver  and  driven  sections.  The  driven 
section  terminates  with  an  end  wall  plate  and  a .077  nm  thick  brass 
diaphragm  to  separate  the  nozzle  from  the  dri.'en  section  prior  to 
performing  an  experiment.  During  operation  of  the  shock  tube,  the 
Incident  shock  wave  reflects  from  the  end  plate,  produces  stagnation 
conditions  for  operation  of  the  nozzle,  and  ruptures  the  secondary 


Fig.  28  SCHEMATIC  OF  THE  12.7  CM  SHOCK  TUBE  GDL  FACILITY 


diaphragm.  Uniform  flow  Is  maintained  through  the  nozzle  for  about  1 
millisecond.  The  nozzle  Is  a minimum  length,  two-dimensional,  contoured 
slit  nozzle  with  a 1 mm  throat  height  and  an  inviscid  core  area  ratio 
of  20.  The  nozzle  width  Is  12.7  cm.  A detailed  drawing  of  the  shock 
tube  end  wall,  nozzle,  and  cavity  is  shown  in  Figure  29. 

The  shock  tube  was  driven  with  cold  helium  in  these  experiments. 

A 5.88%  COg  - 55%  N2  - 39.12%  he  mixture  commercially  obtained  from 

the  Matheson  Company  is  used  as  the  test  gas.  This  gas  mixture  is 

supplied  with  an  accuracy  of  i 2%  of  each  component.  The  total 

impurities,  as  stated  by  Matheson,  are  less  than  200  ppm  and  the  H2 

and  H2O  content  were  each  stated  as  less  than  50  ppm.  Prior  to 

filling  the  shock  tube  with  the  test  gas  mixture,  it  is  evacuated  to 

.3 

less  than  10  torr. 

Experiments  were  performed  for  stagnation  pressures  of  19.47 
to  23.56  atmospheres  and  a range  of  stagnation  temperatures  from  1175 
to  2240  degrees  Kelvin.  The  stagnation  conditions  are  calculated  from 
normal  shock  relationships  assuming  local  thermodynamic  equilibrium 
using  the  method  of  reference  49  and  the  measured  shock  wave  speed. 

Figure  28  shows  the  instrumentation  for  the  measurement  of  the  shock 
wave  speed.  Three  quartz  piezoelectric  pressure  transducers  (manu- 
factured by  PCB  Piezotronics,  Inc.)  and  electronic  counters  are  used 
to  measure  shock  wave  transit  times.  Pressure  is  measured  using  the 
transducer  nearest  the  end  wall. 

Small -signal  gain  measurements  are  made  downstream  of  the 
nozzle  in  the  laser  cavity.  Five  ports  are  spaced  along  the  length 
of  the  cavity  to  enable  longitudinal  GQ  profiles  to  be  made.  GQ  is 
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measured  by  directing  the  beam  of  a low  power  diagnostic  laser  through 
the  cavity  and  measuring  the  increase  in  power  between  no  flow  and  flow 
conditions.  A schematic  of  the  instrumentation  for  measurements  of 
small -signal  gain  is  shown  in  Figure  30.  The  diagnostic  laser  is  a 
Syl vania  Model  948  CO2  laser.  The  beam  is  first  mechanically  chopped 
to  establish  a reference  level  and  then  sampled  (by  means  of  beam 
splitters)  to  determine  wavelength  and  amplitude  stability.  One 
sampled  beam  is  directed  to  a focusing  mirror  which  reflects  it  to  a 
grating  blazed  at  10  micrometers.  This  beam  is  dispersed  by  the 
grating  and  projected  on  to  an  QptEngineering  thermal  imaging  screen 
which  displays  the  individual  line  or  lines  present  and  their  cor- 
responding mode  structures.  The  second  sampled  beam  is  monitored  by 
detector  #1  to  determine  the  amplitude  stability  of  the  diagnostic 
laser.  The  main  probe  beam  passes  through  the  windows  (Irtran  2, 
anti  reflection  coated)  of  the  cavity  and  is  then  diffused  by  reflection 
from  a rough  surfaced  aluminum  flat  to  ensure  coverage  of  the  active 
area  of  the  detector  i! 2.  A narrow-band-pass  filter  is  inserted  in 
the  beam  path  to  prevent  extraneous  radiation  from  invalidating  the 
gain  measurement.  Detectors  #1  and  #2  are  gold-doped  germanium  photo- 
conductive  cells  operated  at  77°K. 

Figure  31  shows  a typical  oscilloscope  trace  of  the  probe  beam 
as  detected  during  an  experiment.  The  lower  trace  is  a 1 OX  amplification 
of  the  upper  trace.  An  upstream  pressure  transducer  triggered  the 
oscilloscope.  The  flat  initial  portion  of  the  oscilloscope  trace  is 
representative  of  a no-flow  condition  in  the  cavity.  The  increase  in 
signal  indicates  the  start  of  flow  in  the  cavity.  (The  upward 
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deflection  is  indicative  of  gain.).  The  rectangular  pulse  is  the 
chopped  signal.  The  reference  is  indicated  by  IQ  and  the  increase  in 
Intensity  by  Al.  The  value  of  GQ  is  given  by  (IQ  + aI)/Iq  ■ exp  (G0L), 
where  L is  the  geometric  path  length  (in  this  case,  12.7  cm  for  a 
single  tranverse  of  the  cavity).  The  gain  measurements  were  made  on 
the  P(28)  and  P(30)  lines.  These  lines  (rather  than  the  usual  P( 20) 
line)  were  used  because  the  diagnostic  laser  was  more  stable  at  these 
transitions  due  to  mirror  misalignment  within  the  laser. 

Typical  pressure  records  near  the  end-wall  are  shown  in  Figure 
32.  The  initial  abrupt  pressure  increase  occurs  when  the  incident 
shock  wave  passes  over  the  pressure  transducer.  The  second  large 
pressure  increase  occurs  as  the  reflected  shock  wave  sweeps  over  the 
gauge.  A pressure  plateau  persists  for  approximately  1 millisecond 
and  is  subsequently  destroyed  by  wave  Interactions.  Two  types  of 
wave  interactions  are  shown  in  Figure  32.  Figure  32a  is  a pressure 
trace  of  conditions  where  the  reflected  head  of  the  expansion  wave, 
generated  by  the  rupture  of  the  primary  diaphragm,  eventually  inter- 
acts with  the  reflected  shock  wave.  The  expansion  wave  lowers  the 
pressure  and  cools  tin  gas  and  thus  terminates  the  test  after  about  a 
millisecond  running  time.  By  comparison.  Figure  32b  represents  the 
interaction  of  the  reflected  shock  wave  with  the  contact  surface  and 
the  subsequent  generation  of  shock  waves  which  are  transmitted  back 
Into  the  stagnation  region.  These  weaker  shock  waves  cause  the 
additional  pressure  and  temperature  increase  in  the  nozzle  reservoir 
region  which  occurs  after  about  1 millisecond  of  useful  test  time. 


Fig.  32  TYPICAL  OSCILLOGRAPH  TRACES  OF  PRESSURE  BEHIND  THE  REFLECTED  SHOCK  WAVE 


Eventually  the  pressure  rapidly  decays,  which  is  indicative  of  the 
arrival  of  the  expansion  wave. 

6.  Arc  Tunnel  Experiments 

Considering  CQ2-N2-H2Q  mixtures,  small -signal  gain  measurements 

have  been  obtained  In  the  HQL  3-Megawatt  Arc  Tunnel  as  a function  of 

distance  from  the  nozzle  throat.  A measured  6Q  profile  is  shown  in 

Fig.  33,  where  it  is  compared  with  the  theoretically  predicted  profile 

obtained  from  the  analysis  of  reference  (8).  Very  good  agreement  is 

obtained.  The  slight  drop  of  the  experimental  data  below  the  theoretical 

curve  at  large  downstream  distances  is  to  be  expected;  the  real  flow 

is  influenced  by  boundary  layer  growth  and  weak  oblique  shock  patterns 

In  the  constant  area  section  whereas  the  theoretical  analysis  assumes 

an  invlscid,  shock-free  flow.  The  existance  of  a weak  shock  pattern 

50  51 

has  been  experimentally  observed  * . In  fact,  a detailed  experimental 

investigation  of  the  fluid  dynamics  of  short,  minimum  length  nozzles 
is  described  in  reference  51. 

Referring  to  F-1;.  33,  the  Initially  rapid  Increase  in  6Q  as 
a function  of  distance  is  due  to  the  rapid  depopulation  of  the  (100) 

C02  level  by  H20.  The  (001)  level  is  also  being  depopulated,  but  at 
a slower  rate.  The  peak  gain  is  reached  when  f^00  essentially  equilibrates 
with  T.  Downstream  of  the  peak,  GQ  decreases  due  to  the  continuous 
deactivation  of  the  (001)  level.  The  different  rates  of  relaxation 
for  the  (001)  and  (100)  levels  are  clearly  evident  in  Fig.  34,  where 
t:.<joretical  results  for  the  two  vibrational  temperatures  and  the 
translational  temperature  are  shown  as  functions  of  distance  along 
the  duct. 
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In  addition,  the  GQ  measured  at  a fixed  station  (1.9  cm  down* 
stream  of  the  nozzle  exit)  Is  shown  In  Fig.  35  as  a function  of  reservoir 
temperature.  Again,  the  measurements  are  compared  with  theoretical 
predictions  based  on  reference  (8)  and  again  reasonable  agreement 
Is  obtained. 

Examining  Fig.  35,  at  lower  temperatures  GQ  Increases  rapidly 
with  Tq  simply  because  the  total  vibrational  energy  In  the  reservoir 
Increases.  However,  the  relaxation  rates  also  Increase  with  temperature. 
Moreover,  for  a nozzle  of  fixed  area  ratio,  the  exit  static  temperature 
increases  as  TQ  Increases,  hence  resulting  In  larger  values  for  N100 
(which  Is  close  to  equilibration  with  T).  For  these  reasons,  GQ 
will  peak  and  begin  to  decrease  If  TQ  Is  made  large  enough.  This 
behavior  Is  clearly  evident  in  Fig.  35.  The  magnitude  and  location 
of  the  peak  depends  on  nozzle  shape  and  size,  pQ  and  mixture  ratio. 

All  of  the  previous  measurements  have  been  made  along 
the  centerline  of  the  flow.  The  minimum  length,  two-dimensional  nozzle 
employed  for  these  measurements  Is  contoured  to  provide  uniform 
flow  at  the  nozzle  exit.  Because  the  flow  field  is  two-dimensional, 
the  fluid  elements  traveling  along  different  streamlines  experience 
different  flow  gradients;  thus  the  nonequilibrium  history  of  each 
streamline  Is  different.  This  prompts  a question  regarding  the 
uniformity  of  GQ  vertically  across  the  nozzle  exit.  In  order  to 
Investigate  this  question,  measurements  of  GQ  were  made  off  the 
centerline  In  a direction  normal  to  the  flow  at  a station  1.9  cm  down- 
stream of  the  nozzle  exit.  These  measurements  are  shown  In  Fig.  36, 
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and  demonstrate  a remarkably  uniform  GQ  distribution  across  the  entire 
Invlscld  flow  region.  These  off-centerline  gain  measurements  < re  be- 
lieved to  be  the  first  of  their  kind,  and  indicate  that  two-dimensional 
effects  are  not  serious  for  such  Internal  nonequilibrium  nozzle  flows. 

7.  Shock  Tunnel  Experiments 

Small -signal  gain  measurements  have  also  been  made  In  C02- 
Ng-He  mixtures  in  the  NOL  12.7  cm  Shock  Tunnel.  The  results  are  shown 
in  Fig.  37,  where  G0  at  a fixed  station  (1.27  cm  downstream  of  the 
nozzle  exit)  Is  plotted  as  a function  of  T . Also  shown  in  Fig.  37  are 
theoretical  predictions  for  the  P(28)  and  P(30)  transitions  of  the 
10.6  u band.  These  are  the  transitions  on  which  the  diagnostic  laser 
was  operating  for  the  present  shock  tunnel  experiments,  as  discussed 
earlier.  Again,  the  theory  based  on  reference  (8)  yields  reasonable 
agreement  with  experiment,  with  the  exception  of  the  higher  temperature 
range  where  the  measured  data  fall  slightly  below  theory.  It  is 
interesting  to  note  that  the  equilibrium  reservoir  conditions  for 
the  present  measurements  Indicate  1 percent  dissociation  for  C02  at 
Tq  * 2000°K,  growing  to  5 percent  at  TQ  * 2250°K.  The  theoretical 
results  shown  in  Fig.  37  assume  a nonreacting,  constant  composition 
mixture  and  hence  do  not  account  for  the  partial  dissociation  of  C02. 
This  contributes  in  part  to  the  over-prediction  of  GQ  at  temperatures 
above  2000°K. 

Small-signal  gain  measurements  for  C02-N2-lle  mixtures  have 

13 

also  been  made  in  a shock  tunnel  by  Christiansen  and  Tsongas  at 

the  University  of  Washington.  Some  of  these  data  are  shown  in  Fig.  38, 
where  GQ  is  plotted  over  a wide  range  of  pQ.  The  measurements  were 
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Flq  37  VARIATION  OF  SMALL-SIGNAL  GAIN  W.iH  RESERVOIR  TEMPERATURE:  COMPARISON  BETWEEN 
THEORY  AND  EXPERIMENT  AT  A STATION  1.27  CM  DOWNSTREAM  OF  THE  NOZZLE  EXIT, 
pe  *=  20  ATM,  X^qj  = 0.0588,  X^e  * 0.3912,  Xf^j  = 0.55 
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rode  at  the  exit  of  a wedge  nozzle  with  throat  height  ■ 2 mm,  total 

* 

Included  angle  =30°,  and  Ae/A  ■ 10.  Also  shown  in  Fig.  38  are 
the  corresponding  theoretical  predictions  from  the  analysis  of 
reference  (8).  Reasonable  agreement  is  obtained  between  theory  and 
experiment. 

8.  Other  Experiments 

A complete  set  of  GOL  Arc  Tunnel  experiments  have  been 
reported  by  Lee  et  al**\  They  give  comparisons  between  theory  and 
experiment  for  small-signal  gain  and  laser  power  for  COp-Ng-He  and 
H20  mixtures.  The  nozzle  arrangement  was  an  array  of  24  Mach  4 nozzles 
of  1 mm  throat  height,  in  contrast  to  the  single  slit  nozzles 
described  in  the  NOL  experiments.  Some  of  their  results  for  small 
signal  gain  are  shown  in  Figures  39-44,  and  some  power  data  are  shown 
in  Figure  45.  In  these  figures, which  are  taken  directly  from  Reference 

O 

10,  the  theoretical  gain  is  obtained  from  the  method  of  Anderson  . 

Fairly  reasonable  agreement  is  obtained.  The  theoretical  curves 

(the  solid  curves)  for  powei  in  Figure  45  were  obtained  from  an 

44 

approximate  analysis  based  on  Rigrod  for  various  values  of  mirror 
absorptivity.  The  actual  mirror  absorptivity  in  these  experiments 
is  not  known,  and  Figure  45  illustrates  the  tenuous  nature  of  power 
extraction  analyses. 

Another  set  of  arc-tunnel  GDL  experiments  were  performed 
52 

by  Howgt 'm,  Roberts  and  Barr  , and  were  compared  with  their  detailed 
numerical  calculations  of  laser  power  output.  Discrepancies  between 
theory  and  experiment  ranged  from  IX  to  100£,  and  again  illustrates 
the  difficulties  involved  in  accurately  predicting  power  extraction 
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Finally,  with  regard  to  conventional  GDI  arc-tunnel  experl- 
53 

ments,  Mcleary  reports  "approximate"  agreement  between  his  theoretical 
calculations,  and  the  arc-tunnel  results  presented  by  Anderson  and 
Winkler12* 

Other  shock  tunnel  experiments  are  reported  by  Klosterman 

p i n 1 A 1 ^ 

and  Hoffman  , Christiansen  and  Tsongas  , Hertzberg  et  al  , Vamos  , 
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Tennant  et  al  , Kuehn  and  Honson  , Blrlukov  et  al  , Sato  and 

Seklguihi55,  and  Gembarzhevskly  et  al5®.  Thoughout  all  of  these, 

comparison  between  theory  and  experiment  ranged  from  poor  to  very  good. 

At  this  point,  the  reader  Is  reminded  that  a conventional 

gasdynamic  laser  operates  with  an  equilibrium  mixture  of  C02,  N2  and 

H20  or  He  in  the  reservoir,  independently  of  the  way  in  which  this 

hot,  high  pressure  mixture  is  obtained.  It  can  be  formed  in  an 

electric  arc  or  by  means  of  shock  heating,  as  previously  described, 

or  it  can  also  be  formed  by  chemical  combustion  as  initially  reported 

by  Gerry1.  Several  combustion-driven  GDI's  have  been  reported  by 

Gerry1,  Tulip  and  Sequin1**,  Yatsiv  et  al1®,  Meinzer2*1,  and  Gabai 
57 

et  al  . In  each  of  these,  comparisons  are  made  with  theory,  and 

again  various  degrees  of  agreement  are  obtained.  Of  particular 

20 

Interest  are  the  extensive  parametric  tests  of  ileinzer  , including 

20 

the  effect  of  H20  on  GDI  gain  all  the  way  to  50%  ll20  content,  ileinzer 

compares  his  experimental  data  with  resLlts  obtained  from  the  analysis 
8 

of  Anderson  , showing  some  reasonable  agreement. 
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A final  remark  Is  In  order  concerning  the  comparison  of  theory 
25  26 

with  experiment.  Anderson  * has  recently  reported  that,  based 
on  the  current  discrepancies  In  the  existing  rate  data,  an  uncertainty 
of  25-30X  can  be  expected  In  any  calculation  o * gasdynamlc  laser 
performance.  These  thoughts  are  given  in  detail  In  Appendix  A of 
the  present  Notes.  Also,  In  addition  to  uncertainties  in  the  kinetic 
rates,'  any  comparison  between  theory  and  experiment  Is  subject  to 
error  due  to  shock  wave  and  viscous  flow  effects  in  the  laster  cavity; 
such  effects  are  generally  not  Included  in  the  existing  theoretical 
analyses. 


E.  INVERSIONS  BEHIND  NORMAL  SHOCK  WAVES 

3 

Basov  and  Oraevskii  suggested  that  population  inversions 
In  molecular  systems  could  be  obtained  by  rapid  heating  of  the  system, 
as  well  as  by  rapid  cooling  as  in  the  case  of  the  conventional  gas- 
dynamic  laser.  One  of  the  most  rapid  ways  of  heating  a gas  is  by 
means  of  a shock  wave.  Therefore,  let  us  examine  the  purely 
vibrational  relaxation  process  behind  a normal  shock  wave  in  C02- 
N2-He  mixtures,  and  assess  what  inversions,  if  any,  occur  in  the 
nonequilibrium  flow.  This  problem  has  been  examined  by  Anderson 
and  Madden  in  References  58-59. 

Consider  the  gasdynamic  flow  model  shown  in  Figure  46. 

We  make  the  standard  assumptions  of  a stationary,  discontinuous  shock 
front  with  frozen  conditions  immediately  behind  the  front,  hence, 
at  this  location  p2,  T2,  and  u2  are  obtained  from  the  standard 
calorically  perfect  gas  equations  for  normal  shocks,  and  ey^  and 
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ev1b  are  equa^  t0  resPect^ve  upstream  values.  In  turn,  these 
quantities  are  the  boundary  conditions  for  the  downstream  nonequilibrium 
flow,  which  Is  solved  by  forward  numerical  integration  of  the  governing 
steady  flow  conservation  equations  as  functions  of  distance  behind 
the  shock  front.  These  equations  are: 

c/f>  ^ ^ ctu  - £ 


Continuity: 


Momentum: 


( *r/f ) £ + + «£=<=> 

Energy:  fi!  T £ +u  **  £ = <=> 


Rate: 


cf  (?vih>s  CCr 

X ~ u- 


a/ it/ 


x 


&4r 

u. 


(41) 

(42) 

(43) 

(44) 


where  a * |xHe  + |{XC(J  + X;J  ),  p ■ysRT,  and  the  notation  Is  standard. 

The  solution  Is  terminated  when  equilibrium  values  of  the  normal  shock 
properties  are  closely  approached;  equilibrium  normal  shock  properties 
for  mixtures  are  known  in  advance  from  Ref.  49.  A detailed 

discussion  of  the  numerical  aspects  of  the  present  study  can  be  found 
In  Ref.  59. 

Numerical  results  for  a typical  case  are  shown  in  Figs.  47  and  48. 

Figure  47  illustrates  temperature  variations  in  the  nonequilibrium  region 

behind  the  shock  front,  and  clearly  shows  the  rapid  equilibration  of  T^b 

with  the  translational  temperature  T,  whereas  In  contrast,  T ..  relaxes 

V1DI! 

more  slowly.  At  a distance  of  4 ran  downstream  of  the  shock  front,  all 
three  temperatures  have  equilibrated  within  one  percent  of  the  final 
equilibrium  temperature,  v/hich  has  been  taken  from  the  results  of  Ref.  59. 
The  results  shown  In  Fig.  47  reflect  a molecular  collislonal  process 
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which  indeed  leads  to  population  Inversions  behind  the  shock  front, 
as  shown  in  Fig.  48.  These  Inversions  occur  between  the  (04°C)  and 
(001)  levels  and  to  a lesser  degree  between  the  (200)  and  (001  levels 
In  COg.  Examining  Fig.  48,  near  the  shock  front  the  inversions  rapidly 
Increase  due  to  the  rapid  population  of  the  excited  levels  of  mode  I 
while  at  the  same  time  the  excitation  of  mode  II  Is  lagging  far  behind. 
However,  the  inversion  soon  peaks  and  begins  to  decrease  farther  down- 
stream as  the  lower  (001)  level  is  substantially  populated.  Note 
that,  for  the  given  upstream  conditions,  the  Inversions  persist  over 
a length  from  1-2  nm  behind  the  shock  front.  The  vibrational  kinetics 
obey  binary  scaling;  hence,  the  spatial  extent  of  the  population 
Inversions  can  be  Increased  or  decreased  by  a propertlonal  decrease 
or  Increase  in  p^  keeping  T^  and  u<|  the  same.  Many  additional 
results  obtained  from  the  present  study  are  discussed  in  Ref.  59. 

How  do  the  laser  properties  of  this  nonequilibrium  shock 
flow  compare  with  those  obtained  by  rapid  expansions?  First,  the 
conventional,  rapid  expansion  gas  dynamic  laser  creates  a population 
Inversion  between  the  (001)  id  (100)  levels  In  CO2  which  subsequently 
lases  at  \ « 10.6  . In  contrast,  the  inversions  shown  in  Fig.  48  between 
the  (04°0)  and  (001)  levels,  and  between  the  (200)  and  (001)  levels, 
would  correspond  to  laser  transitions  at  bO^and  2^,  respectively. 

An  important  parameter  for  gas  lasers  Is  small  signal  gain,  GQ,  defined 
as  dl/I  * GQdz  where  I Is  the  incident  radiation  Intensity  on  a slab 
of  laser  gas  of  thickness  dz,  and  dl  is  the  increase  In  beam  intensity 
after  traversing  the  length  dz.  As  shown  In  Appendix  A of  Ref.  (32), 
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6q  « (X  h ) IN  « (M  /X)  IN,  where  t2j  is  the  spontaneous  radiative 
lifetime  for  a transition  between  the  upper  and  lower  laser  levels, 

M Is  the  corresponding  quantum  mechanical  matrix  element,  and  III  is  the 
population  inversion.  For  CO2,  computed  values  of  M for  the  50m, 

22u,  and  10.6m  transitions  are  in  the  ratio  0.21  x 10“2:0.21  x 10"2;0.34  x 10 
respectively.®0  Also,  the  shock  Induced  population  inversions  shown 
In  Fig.  48  are  approximately  one  order  of  magnitude  smaller  than 
typical  Inversions  created  in  rapid  expansion  through  supersonic 
nozzles.  In  light  of  the  above  numbers,  a comparison  of  G„  ««.  50  u 
and  22y  behind  a shock  wave  with  GQ  at  10.6m  in  a rapid  expansion 

leads  to  (&o)50u/(6o)i0.6Ma  10  4 and  *V22m^  Vi0.6m*  2 x 10  4* 

Clearly,  the  nonequilibrium  region  behind  a normal  shock  wave  in  C02- 

N2-He  mixtures  produces  a low-gain  medium.  A more  detailed  discussion 

and  comparison  of  these  and  other  laser  properties  are  contained 

in  Ref.  59. 

This  study  indicates  that  population  inversions  occur 
behind  a normal  shock  front  due  strictly  to  translation-vibration  and 
vibration-vibration  molecular  energy  exchanges  in  COg-Ng-He  mixtures. 

However,  the  laser  properties  of  this  shock-induced  nonequilibrium 
flow  are  clearly  not  as  promising  as  those  of  gas  dynamic  lasers  operating 
on  the  principle  of  rapid  expansion. 


II.  GASDYNAMIC  LASERS:  STATE -OF-THE  ART 
A.  INTRODUCTION 

High  energy  gas  lasers  are  a modern  reality.  Moreover,  the 
technology  of  these  lasers  is  growing  vigorously  in  at  least  three 
directions,  namely,  gasdynamic,  electric  discharge  and  chemical  lasers. 

The  basic  physical  process  common  to  all  these  lasers  is  the 
competition  between  stimulated  emission  anrl  absorption  of  mono- 
chromatic radiation,  where  the  radiation  (photon)  energy  corresponds 
to  the  difference  between  two  distinct  energy  levels  of  an  atomic  or 
molecular  system.  It  has  been  shown  in  Section  I-B  that  stimulated 
emission  will  win  over  absorption  if  a population  Inversion  exists  between 
the  two  energy  levels,  i.e.,  if  > 0,  and  therefore  laser 

energy  can  in  principle  be  extracted  from  the  gas.  In  electric  discharge 
lasers,  the  uppsr  energy  level  Is  .referential ly  populated  by  collisions 
with  electrons  energized  by  an  electric  field;  in  chemical  lasers, 
the  products  of  highly  energetic  chemical  reactions  are  formed  directly 
In  vibrational ly  or  electronically  excited  states  with  the  upper  levels 
preferentially  populated;  in  gasdynamic  lasers,  an  Initially  hot  gas 
in  thermodynamic  equilibrium  Is  rapidly  expanded  through  a supersonic 
nozzle,  and  an  Inversion  is  formed  by  differential  relaxation  processes 
In  the  nonequilibrum  nozzle  flow.  Of  these  three  types  of  lasers,  the 

1 2 

gasdynamic  laser  has  produced  the  highest  continuous  wave  power  to  date  * . 

The  results  embodied  in  References  1 ,2, N, 8, 9-14, 16-20, 52-57 
represent  (for  the  most  part)  a "first  generation"  of  GDL  technology, 
where  mixtures  of  COg-Ng-HgO  or  He  are  utilized  at  typical  conditions 
of  po«20  atm,  To«1200°K,  h*»l  rnn,  Afi/A*«20,  and  water  mole  fraction 
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Xjj  0*0.01.  (An  exception  Is  the  high  pressure  work  of  Hertzberg^ 

2 13 

and  Christiansen  which  deals  with  pQ  In  hundreds  of  atmospheres.) 

In  light  of  this  previous  work,  the  purpose  of  the  present  section 
Is  twofold:  (1)  To  survey  some  recent  advances  which  contribute  to 

the  "second  generation"  of  GDI  .echnology,  and  (2)  To  present 
experimental  and  theoretical  results  which  typify  theese  advances. 

In  particular,  this  second  generation  of  GDI's  Is  characterized  by 
higher  reservoir  temperatures  and  pressures,  smaller  nozzle  throat 
heights,  larger  area  ratios,  and  larger  HgO  contents.  Moreover, 

In  some  cases  new  gases  and  lasing  transitions  are  being  employed, 
such  as  the  Interesting  work  of  McKenzie®^  with  CO-^-Ap  mixtures. 

As  a final  Introductory  note.  In  the  following  sections  all 
theoretical  values  for  small -signal  gain  and  maximum  available 
energy  are  obtained  from  the  time-dependent  nonequilibrium  nozzle 
analysis  discussed  in  Reference  8,  and  are  calculated  with  the 
computer  code  descr  ied  In  Reference  62.  As  discussed  in  Section  I-C, 
this  analysis  fully  couples  the  vibrational  kinetic  rate  equations 
with  the  governing  quasl-one-dlmensional  flow  equations  of  continuity, 
momentum  and  energy.  Nonequilibrium  conditions  are  Included  both 
upstream  and  downstream  of  the  nozzle  throat.  A simplified  vibrational 
model  Is  used  which  approximates  the  detailed  molecular  energy 
transfer.  The  experimental  results  of  References  9-11  show  that  the 
model  and  analysis  of  Reference  8 yield  reasonable  agreement  with 
experimental  data  for  gasdynamlc  laser  gain.  Also  In  the  present 
discussion,  the  maximum  available  energy,  e , Is  always  the  local 
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value  at  the  nozzle  exit;  hence,  emax,  takes  into  account  the 

vibrational  relaxation  In  the  nozzle.  In  this  manner,  the  present 

e Is  different  from  that  defined  In  Reference  1,  where  the 
max 

mxlmum  available  energy  Is  based  on  reservoir  conditions.  The 
precise  definition  and  calculation  of  eMX  is  given  in  Section  1-C 
of  the  present  notes. 

i 

B.  INCREASED  H20  CONTENT 

For  C02-N2  gasdynamlc  lasers,  the  presence  of  »20  or  He  Is 
extremely  beneficial  for  the  rapid  deactivation  of  the  lower  laser 
level;  indeed,  such  a -catalyst”  is  necessary  for  the  production  of 
high  gain  and  efficient  extraction  of  laser  power.  The  typical 
GDI  described  in  Reference  1 is  combustion  u.-iven,  therefore  «20 
rather  than  He  is  the  meaningful  catalyst  for  such  devices.  Experience 
gained  with  the  first  generation  of  GDL  technology  indicated  that 
gain  is  a very  sensitive  function  of  «20  content;  along  with  the 
beneficial  dead ■ vat ion  of  the  lower  level,  there  is  also  the  competing 
detrimental  deactivation  of  the  upper  level  due  to  «20.  Hence,  with 
the  first  generation  of  GDL's,  on  optimum  amount  of  H20  in  the  gas 
mixture  was  found  to  be  on  the  order  of  1 percent,  i.e..  * 0.01. 

For  H20  contents  much  larger  than  this,  the  collisional  deactivation 
of  the  upper  laser  level  is  overwhelming.  Such  effects  are  shown  by 
the  lower  curves  in  Figures  49  and  50.  In  Figures  49  and  50.  theoretical 
values  for  peak  small-signal  gain  and  maximum  available  energy, 
respectively,  are  given  as  a function  of  X,,^.  (In  the  spirit  of 
the  SI  system  of  units,  eMX  is  given  in  kjoules/kgm;  however,  the 
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Peak  gain  a>  a function  of  H2O  content 
for  first  and  second  generation  GDL's# 


more  connon  GDI  unit  Is  kjoules/lb^,  where  1 kjoules/lbm  - 0.454 
kjoules/kgm.)  The  lower  curves  In  both  these  figures  typify  the 
first  generation  technology,  for  which  the  maximum  peak  gain  is 
clearly  seen  to  occur  around  X,^g  ■ 0.01. 

The  low  value  of  XH  Q*0.01  for  optimum  gain  places  a 

severe  constraint  on  candidate  fuels  for  GDL's;  most  conventional 

fuel -oxidizer  combinations  which  produce  C02  and  N2  also  produce 

H20  far  In  excess  of  1 percent.  (See  Reference  63.)  Therefore,  recent 

efforts  have  been  made  to  study  the  effects  of  higher  H2Q  content 

on  GDL's,  say.  In  the  range  from  1 to  10  percent.  In  this  vein, 

Yatslv  et  al19  and  Tulip  and  Seguln18  have  operated  GDL's  with 

attendant  H20  contents  up  to  8 percent.  (Tulip  and  Seguin  have  very 

recently  used  gasoline  In  a GDL.)  Also,  some  earlier  gain  measurements 

made  In  the  NOL  3-Megawatt  Arc  Tunnel12  covered  a H20  range  from  0 to 

6 percent.  These  experiments  have  shown  that  small -signal  gain  does 

Indeed  persist  at  high  H20  content.  This  is  also  clearly  shown  by 

the  lower  curve  of  Figure  49,  where  GQ  persists  for  H20  content  at 

least  as  high  as  10  percent,  albeit  substantially  reduced  from  its 

optimum  value  which  occurs  around  1 percent.  Of  more  consequence 

for  first  generation  GDL's,  however.  Is  the  drastic  reduction  of 

e shown  In  the  lower  curve  of  Figure  50  for  higher  H20  content, 
max  * 

Therefore,  we  are  led  to  the  conclusion  that,  for  the  h , TQ,  pQ  and 
A /A*  characteristic  of  first  generation  technology,  GDL's  can 
operate  with  high  H20  content,  but  with  severe  penalties  In  efficiency 
I.e.,  the  percentage  loss  In  GQ  and  e^  Is  very  large  as  XH^Q 
Is  Increased  from  0.01  to  0.1. 


94 


On  the  other  hand,  high  values  of  X^  q appear  to  be 

2 £ 

much  more  compatible  with  second  generation  technology  (hower  h , 

A 

higher  TQ,  pQ  and  Afi/A  ).  This  Is  clearly  shown  by  the  upper 

curves  of  Figures  49  and  50,  which  typify  this  n ewer  technology. 

For  A./A*  ■ 50,  h*  ■ 0.3mm,  and  T ■ 1800*K,  the  theoretical  results 
« 0 

shown  In  Figure  49  Imply  thit  peak  gain  Is  maximum  for  XH^p0.025, 
and  that  this  maximum  Is  fairly  flat  for  X^  q from  0.01  to  0.06. 

Indeed,  a more  detailed  study  of  the  theoretical  results  reveals 
that,  for  the  stronger  expansions  associated  with  the  second  gener- 
ation nozzles,  the  tower  laser  level  as  well  as  the  upper  level 
tends  to  freeze  Inside  the  nozzle.  Hence,  more  HgO  Is  necessary 
to  promote  rapid  equilibration  of  the  lower  level  with  the  trans- 
lational energy  of  the  gas.  This  effect  Is  shown  In  Figure  51, 
which  Illustrates  the  variations  of  the  upper  and  lower  laser  levol 
vibrational  temperatures,  TQq^  and  T^qq,  respectively,  as  a 
function  of  distance  for  several  different  values  of  X^q.  Note 
that  for  XH  0 ■ 0.01,  the  lower  level  Is  far  from  equilibrated, 
and  that  values  on  the  order  of  X^Q  * 0.04  are  needed  to  promote  rapid 
equilibration.  Therefore,  we  are  led  ^ conclude  that  second 
generation  GDL's  are  more  amenable  to  higher  HgO  content,  and  that  a 
wider  range  of  fuels  for  GDL's  is  now  a possibility. 

The  gain  profiles  are  markedly  different  between  high  and 
low  H2O  content.  This  Is  shown  In  Figure  52,  where  GQ  is  given  as  a 
function  of  distance  downstream  of  the  nozzle  exit.  In  comparison 
to  the  case  for  X^  q ■ 0.01,  note  in  Figure  52  that  the  peak  GQ  for 
Xu  n * 0.04  occurs  nearer  the  nozzle  exit,  and  that  gain  drops  off 
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more  rapidly  with  distance.  This  Implies  a constraint  on  the  laser 
cavity  design  for  mixtures  with  high  HgO  content,  namely,  that 
power  has  to  be  extracted  In  shorter  distances.  However,  such  rapid 
power  extraction  appears  feasible,  for  example,  with  an  unstable 
resonator,  as  described  by  Gerry\  and  as  discussed  In  a subsequent 
section  of  these  Notes. 

The  results  shown  In  Figures  49-52  are  theoretical,  and  may 

be  a bit  optimistic.  The  kinetic  rates  embodied  In  these  results 

38 

are  obtained  from  Taylor  and  Bltterman  , with  the  lower  laser 
level  deactivation  rate  due  to  HgO  modified  to  account  for  the 
calculations  of  Sharma^.  (See  Section  1-D  of  the  present  Notes.) 

They  do  not  Include  the  up-dated  rates  described  In  References 
25  and  26,  and  reproduced  In  Appendices  A and  B of  the  present 
Notes.  The  effect  of  the  newer  rates  on  such  calculations,  along 
with  a comparison  with  the  experimental  data  of  Vamos^,  is  shown 
In  Figures  53  and  54.  Vanios  has  conducted  a systematic  series  of  shock 
tunnel  tests  covering  1st  generation  (Figure  53)  GDL's  (20:1 
nozzle  area  ratio),  as  well  as  2nd  generation  (Figure  54)  GDL's 
(50:1  nozzle  area  ratio).  In  both  Figures  53  and  54,  the  theoretical 
and  experimental  results  for  small-signal  gain  show  that  GQ  persists 
at  hi gli  H2Q  content.  Therefore,  even  though  the  results  of  Figures 
49-52  may  be  slightly  optimistic,  the  trends  and  conclusions  shown  are 
still  valid. 


C.  INCREASED  TEMPERATURE 

The  obvious  advantage  of  Increasing  the  reservoir  temperature, 
Tq,  of  GDL's  is  that  the  reservoir  vibrational  energy  per  unit  mass 
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Fig.  54 

Variation  of  small- 
signal  gain  with  water  vapor 
content.  Location  1.27  cm 
downstream  of  nozzle  exit. 

Seven  mole  percent  C0-, 
stagnation  conditions*1800*X 
and  37.5  atm.  Water  production 
by  addition  of  U2  and  02. 
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Increases.  On  the  other  hand,  a disadvantage  of  Increasing  TQ 
Is  that  colllslonal  deactivation  becomes  more  rapid,  and 
therefore  not  all  of  the  Increased  reservoir  vibrational  energy 
Is  realized  at  the  nozzle  exit.  Another  disadvantage  Is  Increased 
heating  of  GDL  components.  However,  Increasing  TQ  Is  certainly 
a viable  means  of  Increasing  output  power  per  unit  mass  flow, 

I 

and  high  temperatures  are  therefore  a major  characteristic  of 
second  generation  GDL's.  Because  the  nozzle  exit  temperature 
should  be  low  (Texjt*300*K)  In  order  to  keep  the  lower  laser  level 
population  low,  then  Increasing  TQ  automatically  Implies 
Increasing  A /A*.  For  GDL's  that  exhaust  to  the  atmosphere  with 
a fixed  diffuser  efficiency,  this  also  automatically  Implies 
Increasing  pQ.  Therefore,  In  light  of  the  above  discussion,  the 
reasons  are  clear  for  choosing  the  conditions  which  contrast  the 
first  and  second  generation  technology  In  both  Figures  49  and  50. 

The  purpose  of  this  section  Is  to  discuss  the  results 
of  a numerical  experiment  which  Indicate  the  advantages  of  Increasing 
Tq.  These  results  are  shown  In  Figure  55,  where  GQ  and  emgx  are 
plotted  versus  TQ,  For  each  numerical  point  shown,  Ag/A*  Is 
different;  the  area  ratio  Is  Increased  as  TQ  Increases  In  order 
to  maintain  Texita300°K*  For  the  S0lte  consistency,  the  reservoir 
pressure  Is  held  fixed  at  30  atm  for  these  points,  l.e.,  the  GDL 
Is  assumed  to  exhaust  to  a variable  back  pressure  generally  less 
than  1 atm. 
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Examining  Figure  55,  the  Improvements  brought  about  by 
Increasing  TQ  from  1200  to  1800#K  are  clearly  seen;  peak  G0  Is 
dramatically  enhanced,  and  e^  Increases  by  a factor  of  five. 

This  Is  typical  of  the  Improvements  Incurred  with  second  generation 

technology. 

Now,  let  us  postulate  a third  generation  technology  where 
T further  Increased  to  2400#K  and  Ag/A  *100,  Here,  the 
Improvments  are  not  so  striking;  GQ  Is  beginning  to  plateau  due 
to  Increased  colllslonal  deactivation,  and  the  Increase  In  emax 
Is  less  than  a factor  of  two  between  TQ  * 1800  and  2400%  Also, 
recall  that  all  the  points  shown  In  Figure  55  are  for  pQ  ■ constant  ■ 

30  atm.  If  our  postulated  third  generation  GDL  exhausts  to  1 atm, 
then  pQ  would  have  to  be  higher  than  30  atm  to  overcome  the  Increased 
total  pressure  losses  associated  with  higher  Mach  number  flows.  An 
Increase  In  pQ,  by  Itself,  causes  increased  collisional  deactivation 
In  the  flow,  and  therefore  the  actual  improvements  of  third  generation 
GDL's  will  be  less  than  those  shown  In  Figure  55.  On  the  other 
hand,  for  GDL's  that  exhaust  to  pressures  less  than  1 atm,  the  full 

advantages  shown  in  Figure  55  might  be  realized. 

In  light  of  the  above  discussion,  the  following  conclusions 
are  made.  By  Increasing  TQ  much  above  1800%  only  moderate  improvements 
(factors  of  2 or  3)  are  possible.  However,  these  improvements  can  be 
Important.  Therefore,  a third  generation,  very  high  temperature 
GDL  may  be  worth  further  consideration.  In  the  process,  many  Important 
questions  will  have  to  be  answered;  for  example,  are  there  appropriate 
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fuels  that  yield  a 3000°K  combustion  temperature,  and  how  difficult 

will  It  be  to  maintain  the  critical  elements  of  GDL  hardware  at 

these  temperatures? 

D.  BINARY  SCALING 

1 50 

Gerry  and  Greenberg  et  al  have  pointed  out  that 
current  GDL's  can  be  scaled  according  to  the  parameter  pQh*.  This 
Is  anticipated  because  two-body  collisional  deactivation  dominates 
the  kinetics.  The  purpose  of  this  section  Is  to  examine  In  more 
detail  the  Implications  of  binary  scaling  on  GDL  flows. 

Figures  56  and  57  contain  numerical  results  for  peak  GQ 
and  e(nax>  respectively,  plotted  versus  pQh*  for  a set  of  GDL 
nozzle  flows  with  the  same  mixture  ratio,  nozzle  contour,  T0  and 
A /A*  * 20.  Figures  58  and  59  contain  similar  results,  except 

C 

that  A0/A*  * 50.  In  each  figure,  only  pQ  and  h*  are  varied  as  shown. 

* 

Clearly,  the  peak  GQ  and  emflx  are  unique  functions  of  the  pQh 
product,  thus  explicitly  demonstrating  the  binary  scaling.  In 
addition,  several  trends  are  important  to  note,  as  follows. 

(1)  The  variation  of  GQ  and  e with  pQh*  Is  not 

linear,  l.e.,  halting  the  pQh*  does  not  necessarily  double  the  gain 

* 

or  ematf.  This  Is  particularly  true  for  Gn  at  low  values  of  p h <1  atm-cm, 

inaX  O O 

as  seen  In  Figures  56  and  58. 

(2)  Figures  56  and  58  show  that  P0h*»l  atm-cm  is  a 

reasonable  value  of  the  binary  scaling  parameter  for  high  gain.  This 

50 

value  has  been  stated  earlier  by  Greenberg  . In  the  same  vein, 

Gerry*  has  pointed  out  a criterion  for  adequate  freezing  of  the 
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Fig.  56 

Binary  scaling  of  gain  for  a low 
area  ratio  nozzle* 
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nozzle  flow,  namely,  that  p0h*<2 ( P7^f f ) t*u*  Th1s  criterion  Is  simply 
a statement  that.  In  the  throa t region,  the  flow  residence  time, 
h*/u*,  should  be  less  than  the  effective  relaxation  time,  (Pfeff)j*/ 
p^5(peff)T*/0.5pQ.  For  the  flow  conditions  In  Figures  56-59, 
2(p?£ff)T*u*“l .25  atm-cm;  hence,  a value  of  pQh*»l  atm-cm  satlfles  the 
above  criterion. 

(3)  Major  Increases  In  GQ  can  not  be  obtained  by  reducing 
pQh*  below  1 atm-cm.  This  Is  clearly  seen  In  Figures  56  and  58,  where 
Gq  essentially  reaches  a plateau  for  pQh*<l  atm-cm.  The  basic  reason 
for  this  behavior  is  that,  for  low  values  of  pQh*,  both  the  upper  and 
lower  laser  levels  tend  to  freeze  in  the  nozzle  expansion.  In  fact, 
in  the  theoretical  limit  of  pQ  ■ 0,  hence  pQh*  » 0,  both  the  upper  and 
lower  laser  levels  would  freeze  at  their  equilibrium  reservoir  values, 
no  population  inversion  would  exist  anywhere  in  the  flow,  the  gain 
would  be  negative  and  e would  be  zero.  These  trends  are  illustrated 

II  IQ  A 

In  Figure  60,  which  gives  the  vibrational  and  translational  temperatures 
for  pQh*  * 0.21  atm-cm  (solid  lines)  and  1.2  atm-cm  (dashed  lines) 
as  a function  of  distance  from  the  nozzle  exit.  First,  consider  the 
solid  lines.  For  pQh*  « 0.21  atm-cm,  the  tendency  for  T1Q0  to 
"hang  up"  is  quite  evident.  Indeed,  a downstream  distance  of  more  than 
30  cm  Is  required  before  T-jqq  equilibrates  with  T.  At  30  cm,  the  peak 
gain  of  1.28  m"1  is  finally  achieved.  In  comparison,  consider  the 
dotted  lines  for  pQh*  « 1.2  atm-cm.  The  upper  level  vibrational 
temperature,  Tqq1,  is  less  than  the  previous  case;  by  itself,  this 
will  reducf  GQ.  However,  Tlfl0  has  relaxed  quite  rapidly  and  is 
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essentially  equilibrated  near  the  nozzle  exit;  this  will  favor  the 
gain.  For  this  case,  a peak  GQ  of  1.2  is  achieved  near  the  nozzle 
exit. 

With  these  trends  In  mind,  there  appears  to  be  little  benefit 
in  striving  for  exceptionally  small  values  of  pQh*  In  gasdynamlc 
lasers.  Indeed,  It  appears  to  be  sufficient  to  satisfy  the  criterion, 

P0h*<2tpVf)T*u*- 
E.  CARBON  MONOXIDE  GPL 

In  our  discussion  of  modern  advance  In  gasdynamlc  lasers,  we 
would  be  remiss  not  to  mention  specifically  the  pioneering  work  of  McKenzie 
with  carbon  monoxide.  All  previous  sections  of  these  Notes  have  dealt 
exclusively  with  COg-.^-HgO  or  Cl^-i^-He  mixtures,  with  the  lasing 
transition  at  10. 6u.  However,  this  Is  by  no  means  the  only  possibility 
for  gasdynamlc  lasers.  McKenzie  has  reported  independent  experimental 
and  theoretical  results  for  population  Inversions  and  laser  power 
extraction  Irt'rapldly  expanding  CO-t^Ar*  mixtures.  The  details  are 
clearly  explained  In  Reference  61,  and  hence  no  elaboration  will  be 
given  here.  However,  the  following  Important  points  are  emphasized. 

The  molecular  kinetics  of  CO-Ng-Ar  mixtures  rely  upon  the  freezing 
of  vibrational  energy  In  and  CO,  both  of  which  have  unusually 
long  relaxation  times.  Furthermore,  vibrational  energy  is  pumped 
from  NgC v»l ) to  C0(v»l),  similar  to  the  conventional  COg-?^  system. 

However,  the  similarity  stops  there,  because  tne  vibrational  energy 
Is  rapidly  distributed  upward  to  the  higher  lying  levels  of  CO  due 
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to  rapid  vibration-vibration  (V-V)  transfer.  The  anharmonlclty 

of  these  upper  levels  promotes  the  rapid  V-V  transfers,  as  discussed 
64 

by  Treanor  et  al  . The  net  result  Is  a non-Boltzmann  distribution 
of  molecules  among  the  upper  CO  vibrational  levels.  This  distribution 
may  lead  to  a vibrational  population  Inversion  In  the  usual  sense  (more 
particles  In  an  upper  vibrational  level  than  a lower  level);  or 
more  likely  It  will  result  In  a partial  Inversion  for  a given  rotation- 
vibration  transition,  even  though  a total  vibrational  Inversion  does 
not  exist.  Moreover,  the  CO  gasdynamic  laser  generally  emits 
simultaneously  at  a number  of  wavelengths  In  the  5p  range.  All  these 
points  may  add  up  to  an  order  of  magnitude  Improvement  of  CO  over 
conventional  C02-N2  GX  systems  In  terms  of  power  output  per  unit 
enthalpy  flux  of  the  flow61. 

The  CO  gasdynamic  laser  suffers  from  one  practical 
disadvantage  — the  requirement  for  high  Macl.  number  flows.  The 
population  inversions  In  CO  are  enhanced  by  very  low  translational 
temperatures  of  the  gas,  on  the  order  of  60°K.  In  order  to  achieve 
such  a low  temperature,  the  flow  must  be  expanded  In  supersonic  nozzles 
with  area  ratios  (exit  to  throat  area)  on  the  order  of  1000  — a 
factor  of  20  larger  than  conventional  GDI's.  To  an  aerodynamlclst, 
this  Implies  high  Mach  number  flows  with  large  total  pressure  losses. 
Since  many  practical  applications  of  GDL's  Involve  exhausting  the  gas 
to  the  atmosphere,  then  either  extremely  high  reservoir  pressures  are 
required,  or  a dramatic  breakthrough  In  supersonic  diffuser  technology 
Is  necessary.  This  amounts  to  a serious  limitation  on  the  practical 
use  of  such  lasers. 


F.  DOWNSTREAM  MIXING  6ASDVNAHIC  LASERS 


How  do  we  squeeze  the  most  out  of  the  g,  dynamic  laser 
concept?  The  answer  to  this  question  may  rest  In  some  recent  work  on 
downstream  mixing  GDL's,  which  represents  the  frontier  of  gasdynamlc 
laser  research  today. 

Consider  a conventional  C02-N2-H20  or  He  GDI  as  discussed  In 
these  Notes.  Emphasis  Is  made  that  the  primary  reservoir  of  energy 
for  laser  extraction  Is  In  the  vlbratlonally  excited  Ng,  even  though 
the  laser  emission  Itself  occurs  via  the  CO2  molecules.  Of  course, 
one  purpose  of  the  supersonic  nozzle  expansion  Is  to  freeze  as  much  of 
this  vibrational  energy  as  possible.  Pure  N2  has  an  extremely  long 
relaxation  time;  hence  It  readily  freezes  In  a rapid  expansion.  However, 
In  a conventional  GDL  where  both  C02  and  H20  (or  He)  are  present  In 
addition  to  N2  In  the  nozzle  expans Ion, such  molecules  act  as  a contaminant 
to  the  U2.  As  a result,  the  C02-N2-H2Q  or  He  mixture  relaxes  faster 
than  pure  U2,  giving  lower  vibrational  temperatures,  hence  lower 
available  laser  power.  This  Is  shown  schematically  In  Figure  61.  In 
addition  to  this  consideration,  to  obtain  maximum  laser  power,  the 
reservoir  temperature,  TQ,  should  be  made  as  high  as  possible.  The 
higher  the  temperature,  the  more  vibrational  energy  Is  present  In  the 
gas  and  subsequently  more  laser  power  can  be  obtained.  However,  In 
conventional  GDL's,  TQ  Is  limited  to  values  around  2300°K;  above  this 
temperature,  C02  begins  to  chemically  dissociate,  and  laser  performance 
begins  to  deteriorate.  This  is  In  contrast  to  the  dissociation  of  !12, 
which  occurs  above  4000#K. 
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In  light  of  the  above,  the  rader  can  clearly  see  two  major 
advantages  In  constructing  a gasdynamlc  laser  which  expands  pure 
Ng  through  a supersonic  nozzle,  and  then  mixes  COg  and  HgO  or  He 
downstream  of  the  nozzle  exit.  Here: 

(1)  Reservoir  temperatures  can  exceed  4000°K,  and 

(2)  More  efficient  freezing  of  the  Ng  vibrational  energy  Is 
obtained. 

Such  a device  Is  called  a downstream  mixing  gasdynamlc  laser,  and  Is 
shown  schematically  In  Figure  62.  If  the  Ideal  situation  of  Instantaneous 
mixing  Is  assumed,  then  for  the  conditions  shown  In  Figure  62,  a non- 
equilibrium calculation  using  the  method  of  Reference  62  Indicates  a 
maximum  available  power  of  over  300  KJ/1bm.  This  Is  an  order  of  magnitude 
Improvement  over  the  existing  2nd  generation  GDL's  discussed  In  Section 
II  above.  In  fact,  this  number  Is  compared  with  other  laser  concepts 
In  Figure  63,  which  clearly  Indicates  why  there  is  current  Interest  In 
such  a device. 

However,  there  are  obvious  problems.  Mixing  does  not  occur 
Instantaneously,  and  throughout  the  finite  mixing  region  of  a real  flow, 
vibrational  deactivation  will  constantly  decrease  the  available  power. 

The  question  Is:  Does  enough  deactivation  occur  during  mixing  to  negate 

the  whole  concept?  Work  Is  being  done  to  answer  this  question,  as 
described  below.  There  are  other  problems.  For  example,  can  the 
nozzles  be  made  of  materials  that  can  withstand  such  high  temperatures 
in  the  small  nozzle  throats?  How  can  we  produce  \\^  at  4000°K  on  a 
continued,  practical  basis?  Also,  recall  that  high  Mach  number  nozzles 


Fig.  62 


#r*  required  to  Insure  that  the  gas  translational  temperature 
(hence  the  lower  laser  level  population)  In  the  cavity  will  be 
low,  Hence,  can  diffusers  be  developed  that  will  allow  exhausting  the 
flow  to  the  atmosphere?  Can  the  flows  be  mixed  without  Inducing 
strong  shock  waves  or  major  turbulence  In  the  laser  cavity  that  degrade 
beam  quality? 

The  advantage  of  high  power  output  obtained  In  a downstream 
mixing  GDL  encourages  work  on  this  concept,  even  In  the  face  of  the 
above  problems.  Some  early  experimental  work  has  been  performed  by 

cc 

Bronfln  et  al  using  arc-heated  Ng  with  COg-He  Injected  vertically 

Into  the  supersonic  flow  In  the  cavity.  Their  results  were  the 

first  to  demonstrate  the  potential  of  the  downstream  mixing  concept. 

However,  this  work  was  not  pursued,  presumably  due  to  concern  about 

shock  waves,  hence  poor  beam  quality,  In  the  cavity,  .lore  recently, 

downstream  mixing  GDL  experiments  using  tangential  mixing  have  been 

reported  by  Hilewskl  et  al®®  in  a shock  tube.  (Note:  Reference  68 

52 

contains  more  recent  data  than  Reference  67.)  Also,  Howgate  et  al 
report  some  similar  measurements  with  COg  Injected  both  at  the  nozzle 
exit  and  upstream  of  the  throat.  Nov/ever,  the  most  complete 
Investigation  by  far  has  been  conducted  by  the  team  of  Borghl  and 
Charpenel®0,  who  have  performed  detailed  vibrational  nonequillbrium 
flow  field  calculations  of  tangential  downstream  mixing  using  a finite- 
difference  solution  of  the  turbulent  boundary  layer  equations,  and  of 
Taran,  Charpenal  and  Borghl*0,  who  report  a systematic  series  of  arc- 
tunnel  experiments  complementing  the  theoretical  predictions.  Some 
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of  the  results  of  Taran  et  al70  are  shown  In  Figures  64-68.  They 
examined  mixing  downstream  of  the  nozzle  exit  (Figure  64)  as  well 
as  mixing  In  the  throat  region  (Figure  65).  Theoretical  results  for 
9a1n  profiles  are  shown  In  Figure  66  for  mixing  downstream  of  the 
exit,  and  In  Figure  67  for  mixing  In  the  throat  region.  Clearly, 
mixing  In  the  throat  region  shows  higher  gains  by  an  order  of 
magnitude.  This  should  not  be  construed  as  a general  result;  rather. 
It  IS  due  to  the  investigator's  choice  of  mixing  a high  Mach  number 
stream  of  N2  with  a low  Mach  number  stream  of  C02,  as  shown  In 
Figure  64.  As  a result,  a considerable  amount  of  the  high  total 
temperature  of  the  «2  stream  was  recovered  In  the  cavity,  leading  to 
high  static  temperatures  (800-1000'K) , hence  low  gains.  This  choice 
of  mixing  high  and  low  Mach  number  streams  was  made  to  enhance  the 
mixing  process.  Such  a high  temperature  recovery  was  not  obtained 
in  the  scheme  shown  in  Fig.-re  65.  These  theoretical  results  were 
reasonably  confirmed  by  arc-tunnel  results,  exemplified  by  the 
experimental  data  shown  In  Figure  68.  The  reader  Is  urged  to  study 
References  69  and  70  for  a full  exposition  of  these  studies.  They 
represent  the  frontier  of  the  state-of-the-art.  In  fact,  for  the 
reader's  convenience.  Reference  70  Is  reproduced  as  Appendix  C of  the 
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Fig.  67  - Gala  distribution  on  th*  P 16  11m  for 
nixing  at  tba  throat.  Tba  stagnation  conditions 
for  »2  arc  2500*1  and  15  At*.  Tba  J#  COj,  7#  Ha 

■lxtura  la  at  rooa  taeperatura  t tba  half  slot 
Might  of  tba  Injector  la  0.3  am  The  throat 
Might  for  tha  I?  (between  tba  injector  and  ooa 

notsla  wall  la  0.7  ej)»  now  pro  portion!  ora 
4*  S„  Y#  Ha.  IT*  CO,. 


HI  OPTICAL  CAVITY  CONSIDERATIONS  AND 
FLUID  DYNAMIC  PROBLEMS 


A.  INTRODUCTION 

To  this  point  we  have  discussed  the  many  facets  of  GOL's 
which  bear  on  population  inversion,  small-signal  gain  and  maxi™  available 
power.  To  a lesser  extent,  some  of  the  aspects  of  power  extraction  were 
considered  in  Section  IC.  However,  power  extraction  is  more  than  Just 
solving  the  radiative  transfer  equation.  dl/I  = t dy  between  mirrors  in 
the  cavity,  coupled  with  the  nonequilibrium  supersonic  flow.  It  is  also 
a matter  of  considering  the  optical,  design  of  the  cavity,  taking  into 
account  the  actual  path  of  the  light  rays  (geometric  optics),  mirror 
losses  (such  as  absorption),  diffraction  (physical  optics),  and  the 
density  Inhomogeneities  of  the  flowing  gas  which  cause  phase  distortions 
in  the  laser  beam.  In  essence,  the  basic  philosophy  of  optical  cavity 
design  is  stated  by  Clark71:  "Achieve  the  maximum  output  power  in  a 

beam  with  minimum  divergence,  limited  only  by  diffraction."  Such 

considerations  are  the  subject  of  this  section. 

Also,  emphasis  is  made  that  the  optical  phase  distribution 

across  the  output  beam  from  a laser  cavity  should  be  constant  (uniphase). 
and  that  the  beam  divergence  angle  should  be  as  small  as  Nature  allows 
(limited  only  by  diffraction  through  the  output  aperture).  For  these 
conditions,  the  laser  beam  can  be  focused  to  a spot  of  minimum  diameter. 
Any  deviation  from  these  conditions  will  reduce  the  beam  intensity  in 

the  far  field. 
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8.  CAVITY  TYPES 
1.  Amplifiers 

The  first  generation  of  GDL's  have  extracted  power  by  using 
the  cavity  as  an  amplifier,  such  as  shown  In  Figure  69.  Here,  an 
independent  beam  from  a CC^  laser  (usually  a “conventional"  electric 
discharge  laser)  Is  fed  Into  the  GDL  cavity  through  an  Input  aperture. 

This  CO2  laser  Is  called  the  Master  Oscillator  (MO).  The  beam  Is  then 
reflected  back  and  forth  Inside  the  GDL  cavity  for  several  passes, 
being  amplified  by  the  laser  medium  on  each  pass.  On  the  last  pass, 
the  beam  leaves  the  cavity  either  through  small  holes  In  the  output 
mirror  or  through  some  type  of  window.  At  present,  extracting  the 
last  pass  through  a solid  window  Is  impractical  because  current  window 
materials  can  not  tolerate  the  high  heating  loads  Imposed  by  the  high 
pov/er  beam.  This  leads  to  the  use  of  aerodynamic  windows,  to  be  discussed 
later.  Referring  again  to  Figure  69,  the  cavity  Is  acting  as  a power 
amplifier  (PA),  and  the  whole  arrangement  of  oscillator  and  cavity  Is 
called  the  Master  Oscillator  Power  Amplifier  (MOPA).  The  number  and 
arrangement  of  passes  through  the  cavity  depend  on  the  designer’s 
intentions;  several  possible  geometries  are  shown  in  Figure  69,  which 
Is  taken  from  Reference  71. 

There  are  three  major  advantages  of  the  MOPA  for  extracting 

power: 

(a)  The  Imput  beam  from  the  Master  Osclllat;'",  assuming 
It  Is  properly  designed,  Is  already  in  constant  phase.  Hence,  If  the 


Typical  amplifier  folding 
geometries . 
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amplifier  medium  In  the  cavity  Is  homogeneous,  then  the  output  beam 
should  be  of  constant  phase  also. 

(b)  The  power  level  of  the  output  beam  of  a 110 PA  can  be 
readily  Increased  or  decreased  simply  by  controlling  the  Input  beam 
from  the  MO.  This  adds  a certain  flexibility  to  GDL  operation. 

(c)  The  timewise  (temporal)  variation  of  the  output  beam  can 
also  be  controlled  In  a similar  manner. 

Referring  again  to  Figure  69,  a disadvantage  of  the  MOPA  Is 
that  most  of  the  power  available  In  the  supersonic  gas  In  the  cavity 
Is  extracted  only  during  the  last  beam  pass;  the  role  of  the  other  passes 
Is  simply  to  build  up  the  Intensity  so  that  the  last  pass  reaches 
saturation  Intensity.  At  saturation  Intensity,,  the  beam  Is  extracting 
energy  from  the  C02  molecules  just  as  fast  as  the  il2  can  pump  energy 
Into  the  upper  laser  level  (001)  of  C0?.  In  the  meantime,  during  the 
earlier  beam  passes  that  occur  upstream  of  the  last,  saturating  pass, 
available  laser  energy  in  the  gas  Is  being  lost  due  to  colllsional 
deactivation  of  the  molecules.  Hence,  a IKJPA  arrangement  Is  in  general 
an  Inefficient  way  of  extracting  the  available  laser  energy. 

2.  Resonators 

A resonator  Is  a self-contained  extraction  device,  where  the 
laser  radiation  is  Initially  triggered  by  random  spontaneous  emission 
and  then  rapidly  Increased  by  stimulated  emission  In  a direction 
perpendicular  to  mirrors  on  both  sides  of  the  cavity.  Three  types  of 
resonators  are  shown  In  Figure  70,  taken  from  Reference  71.  Resonators 
have  the  advantage  of  rapidly  building  up  saturation  Intensity  within 
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the  cavity,  hence  being  more  efficient  In  general  than  the  MQPA 
discussed  above.  They  have  the  disadvantage  of  allowing  no  Intensity 
control,  and  the  optical  alignment  is  more  difficult  than  IKJPA's. 

However,  their  efficient  operation,  which  Is  related  to  their  ability 
to  extract  energy  In  a short  flow  length,  makes  such  resonators  of 
extreme  Interest  In  GDI  design. 

There  are  basically  three  types  of  resonators,  and  they  all 
have  been  used  at  one  time  or  another  on  a GDL.  These  types  are 
shown  In  Figure  70.  The  major  distinction  between  them  Is  that 
light  rays  which  are  slightly  off-axis  tend  to  "walk  out"  of  the 
plane-parallel  and  unstable  resonators,  but  are  contained  by  the 
stable  resonator. 

Diffraction  plays  a very  important  role  in  determining  what 
type  of  intensity  pattern  exists  in  the  output  laser  beam  which  emerges 
from  the  output  mirror.  As  stated  by  Clark^,  to  produce  a phase 
coherent  beam,  all  parts  of  the  laser  cavity  should  share  energy  through 
diffraction.  In  this  way,  the  optical  phases  of  all  the  emitting  molecul 
are  locked  together  through  cross-coupling  within  the  cavity,  and  a 
laser  beam  of  the  lowest  order  transverse  mode  will  be  produced.  Some 
cavities  do  this  job  better  than  others.  For  example,  a number  of 
laser  transverse  modes  can  oscillate  within  a cavity,  each  producing 
a different  intensity  pattern  in  the  output  beam.  The  Intensity 
patterns  associated  with  twelve  of  these  modes  are  shov/n  in  Figure  71, 
which  Is  taken  from  Reference  72.  The  best  quality  laser  beam  is  one 
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when  a single  node  is  operating,  with  that  node  being  the  lowest  order 
TEH  node.  (TEII  means  transverse  electromagnetic  mode.)  This  rode 
produces  the  single  spot  seen  in  the  upper  left  comer  of  Figure  71. 
Unfortunately,  it  is  very  difficult  to  obtain  single  node  operation; 
Indeed,  nost  GDI’s  have  in  the  past  operated  multi-mode. 

Mhat  is  the  problem  that  causes  multi-mode  operation?  First, 
consider  a stable  resonator,  and  define  the  Fresnel  number.  Nf.  as 
u . r2/dXi  where  r • mirror  radius,  d • distance  between  mirrors  and 
X is  the  laser  wavelength.  Stable  resonators  with  high  and  low  Fresnel 
numbers  are  shown  in  Figure  72.  Diffraction  effects  in  the  long  cavity 
(the  low  Fresnel  number  cavity)  are  stronger;  this  can  be  seen  by 
inspection  of  Figure  72  where  a slightly  off-axis  ray  effectively  has 
a long  enough  distance  between  mirrors  to  cross-couple  with  large  regions 

of  the  cavity. 
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Hence,  a low  Fresnel  number  cavity  has  a good  chance  of  operating 
with  the  lowest  order  mode.  On  the  other  hand,  again  looking  at 
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Figure  72,  the  opposite  Is  true  of  a high  Fresnel  number  cavity. 

This  Is  also  seen  In  Figure  73,  taken  from  Reference  71.  For  a large 
Np,  the  stable  resonator  can  set  up  several  Independent  modes  of 
operation,  with  little  cross -coupling  by  diffraction.  Unfortunately, 
to  handle  the  large  mass  flows  characteristic  of  GDL's,  the  cavity 
has  to  be  a fairly  high  Fresnel  number  cavity.  A partial  solution 
Is  to  use  an  unstable  resonator  as  shown  In  Figure  73.  Even  though 
Np  may  be  large,  the  rays  of  an  unstable  resonator  spread  out 
within  the  cavity,  and  provide  the  desirable  strong  cross -coupling. 
Hence,  In  an  unstable  resonator,  It  Is  more  difficult  for  the 
higher  order  transverse  modes  to  oscillate.  For  more  details,  the 
reader  Is  urged  to  consult  Reference  71. 

B.  CAUSES  OF  SEA11  DISTORTION  IH  GDL's 

All  of  the  above  discussion  has  assumed  that  the  laser  medium 
(the  supersonic  flow  In  the  cavity)  is  homogeneous.  If  the  medium  is 
not  homogeneous,  then  additional  sources  of  phase  distortion  are 
present.  These  additional  sources  are  examined  in  the  present  section. 

1 . Shock  Waves 

There  are  two  principal  sources  of  shock  waves  in  the  GDL 
cavity,  and  these  are  shown  in  Figure  74,  taken  from  Reference  73. 

The  nozzle  of  a GDL  is  really  a bank  of  small  nozzles  in  a row  — 
small  nozzles  to  promote  vibrational  freezing  in  the  expansion.  This 
Is  shown  in  the  top  sketch  of  Figure  74,  which  looks  down  on  the  top 
view  of  a GDL.  Here  we  see  the  minimum  length  contour  of  such  nozzles, 
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as  well  as  a diamond  shock  pattern  coming  from  the  nozzle  tips, 
now  is  from  left  to  right,  and  the  laser  beam  Is  transverse  to  the  flow 
as  shown.  The  flow  at  the  exit  of  each  nozzle  Is  not  exactly  parallel; 
this  Is  due  to  truncation  of  the  nozzles  to  obtain  structural  strength 
at  the  tips,  and  due  to  the  boundary  layer  growth  along  the  nozzle  walls 
as  well  as  the  wake  growth  downstream  of  the  tips.  Hence,  a diamond 
shock  pattern  Is  always  a characteristic  of  such  multi-nozzles, 
and  was  observed  long  before  the  advent  of  gasdynamlc  lasers  . A 
detailed  survey  of  these  disturbances,  as  well  as  other  fluid  dynamic 
aspects  of  gasdynamlc  lasers,  can  be  found  In  the  recent,  excellent 

survey  by  Russell  in  Reference  75. 

This  diamond  shock  pattern  is  generally  weak.  Moreover, 

If  the  laser  beam  cuts  across  these  shocks  In  the  direction  shown 
in  the  top  of  Figure  74,  then  the  optical  path  will  essentially  be 
the  same  for  each  ray  and  phase  distortion  will  be  minimized.  Hence. 
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the  shocks  emanating  from  the  nozzle  tips  are  not  a major  problem. 

This  Is  not  the  case  for  the  shocks  that  come  from  the  top 
and  bottom  walls  of  the  nozzle  array.  Such  shocks  are  shown  In  the 
bottom  sketch  of  Figure  74.  Here  we  are  looking  at  a side  view  of  the 
GDL,  where  the  nozzle  entrance,  throat  and  exit  are  shown  as  the 
vertlca1  lines.  Flow  Is  still  from  left  to  right.  Consider  the 
expansion  In  the  nozzles.  The  pressure,  density  and  temperature  drop 
so  rapidly  In  the  nozzle  throat  region  that  the  boundary  layer  displacement 
thickness  on  the  flat  top  and  bottom  walls  experiences  an  almost  discontinuous 
Increase.  This  appears  as  a compression  comer  to  the  flow,  and  triggers 
the  two  oblique  shocks  emanating  from  the  top  and  bottom  walls  shown 
In  Figure  74.  Here,  the  rays  of  the  laser  beam  are  parallel  to  the 
plane  of  the  shocks,  and  the  optical  length  of  the  rays  can  vary 
substantially  across  the  beam,  depending  on  the  density  of  the  different 
regions  of  the  shock  pattern.  As  a result,  these  top  and  bottom 
shocks  can  cause  substantial  phase  distortion. 

There  are  two  solutions  to  this  problem.  One  Is  gasdynamic, 

73 

and  Is  discussed  by  Simons  . It  Involves  the  simple  aspect  of 
contouring  the  top  and  bottom  walls  In  the  nozzle  throat  region  In 
such  a fashion  as  to  cancel  the  displacement  thickness  growth.  This 
Is  Illustrated  In  Figure  75,  taken  from  Reference  73.  However,  precise 
contouring  of  the  top  and  bottom  walls  is  a serious  fabrication 
complication  for  the  construction  of  GDL's.  A compromise  is  discussed 
by  Director^,  v/iio  shows  experimentally  that  the  oblique  shock  waves 
can  be  sufficiently  mitigated  by  simply  diverging  the  top  and  bottom 
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walls  by  a constant  small  angle  on  the  order  of  a degree  or  so. 

The  other  solution  Is  optical.  Hoffman  and  Jones77  have  suggested 
the  phase  canclllatlon  optical  configuration  shown  In  Figure  76. 


Fig.  76  Ph«»e  Cancellation  Optical  Configuration 


This  is  characterized  by  a symsetrlc  beam  path  arrangement  wherein  each 
ray  cuts  the  shock  waves  an  equal  number  of  tines  and  each  pass  traverses 
the  flow  at  the  same  angle.  St  an  arrangement  applies,  of  course, 
to  the  HOPA  type  of  optical  cavity  design;  It  Is  not  really  applicable 
to  resonators.  This  phase  cancellation  Idea  has  led  to  an  order  of 
magnitude  Improvement  In  beam  quality77,  and  Is  Illustrated  In  Figure  77, 

obtained  from  Reference  75. 

As  a result  of  both  the  above  solutions,  the  serious  beam 
quality  problems  originally  encountered  with  the  cavity  shock  waves  have 
gradually  been  minimized.  In  fact.  In  a properly  designed  "clean  flow" 
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nozzle  and  cavity,  random  turbulence  In  the  flow  may  be  the  limiting 
factor  on  beam  quality  instead  of  the  shock  waves. 


Fia  77  F,*k  <,r  <i,ld  lBt,n,ltF  •* 

•'  10. 6u  calculated  fro*  0.7u 

IntcrfcroRrama  (Iron  Rafar- 
anca  31). 

2.  Turbulence 

There  is  not  much  that  can  be  said  about  turbulence  at  the 
present  time,  and  its  quantitative  effect  on  phase  distortion.  Given 
the  extent  and  spectrum  of  the  turbulence  of  the  flow  in  the  cavity, 
the  analysis  of  Sutton78  allows  a calculation  of  its  effect  on  the 
beam  in  the  far  field.  However,  a prediction  of  the  turbulence  level 
generated  withir,  the  boundary  layers  and  wakes,  and  that  which  propagates 
downstream  from  the  reservoir  and  subsonic  po. tions  of  the  nozzles, 

Is  difficult  tu  obtain.  However,  we  can  catagorically  state  that 
random  turbulence  will  always  create  a temporal  phase  change  distortion 
in  the  beam,  and  that  such  turbulence  can  be  the  fundamental  limiting 
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factor  In  obtaining  the  best  beam  quality  . 

3.  Non-uniform  Gain  Profiles 

Recall  that  small -signal  gain  in  a GDL  varies  with  distance 
as  shown  In  Figure  33,  and  as  sketched  in  Figure  78  for  cases  with  and 
without  power  extraction.  Clearly,  gain  varies  through  the  cavity  and 
across  the  laser  beam.  Such  a non-uniform  gain  profile  ultimately 
leads  to  a non-uniform  Intensity  variation  across  the  output  beam. 

This  Is  shown  In  Figure  79,  obtained  from  Reference  71.  Here,  the 
Intensity  distribution  of  the  Input  beam  In  a M3PA  (a  Gaussian  profile) 

Is  compared  with  the  corresponding  output  intensity  distribution.  Note 
that  the  upstream  edge  of  the  output  beam  (where  gain  Is  high)  peaks 
to  a very  high  value.  These  results  clearly  show  the  non-uniformity 
of  the  output  beam  as  caused  by  a non-uniform  gain  profile. 

The  results  shown  in  Figure  79  were  computed  based  on  geometric 
optics.  However,  in  the  region  of  peal;  intensity  on  the  upstream 
side,  diffraction  is  important.  Figure  80,  also  obtained  from  Reference  71, 
shows  the  same  case,  except  a more  complicated,  physical  optics 
calculation  is  used  to  obtain  the  output  intensity  distribution.  The 
effect  of  diffraction  is  to  smooth  out  the  sharpness  of  the  peak  as 
predicted  from  geometric  optics  in  Figure  79.  However,  Figure  80  stil’i 
shows  a very  non-uniform  intensity  distribution  across  the  beam. 

The  amount  of  power  contained  in  a focused  beam  in  the  far 
field  does  not  seem  to  be  materially  reduced  by  a non-uniform  distribution 
In  the  near  field.  However,  the  peak  levels  of  the  non-uniform  intensity 
distribution  may  cause  local  overheating  of  the  cavity  output  mirror 
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Geometrical  optics  cal- 
culation of  gain  saturation 
induced  distortion  in  a 
GDL  amplifier. 
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or  window,  which  could  cause  major  problems.  Therefore,  such  non- 
uniform  peaks  as  shown  In  Figure  80  must  be  accounted  for  In  mirror 
and  window  design. 

C.  AERODYNAMIC  WINDOWS 

The  preceding  paragraph  suggests  that  problems  may  arise  with 
mirrors  and  windows  In  high  power  lasers.  In  fact,  how  do  we  physically 
extract  the  high  Intensity  laser  beam  from  the  cavity,  and  allow  It  to 
pass  to  the  surrounding  atmosphere?  One  technique  was  mentioned  In 
Sections  IC  and  D,  namely,  to  have  a number  of  holes  In  the  output  mirror 
through  which  a portion  of  the  radiation  can  pass.  This  Is  effective 
In  extracting  raw  power,  but  causes  major  problems  with  beam  quality. 
Incidentally,  air  leaks  from  the  outside  Into  the  cavity  through  these 
holes  can  have  a deleterious  effect  on  the  supersonic  flow  In  the  cavity 
and  diffuser.  An  alternative  technique  Is  to  simply  extract  the  beam 
through  a solid  window.  However,  a certain  amount  of  the  laser  energy 
Is  always  absorbed  by  the  window  material,  and  for  high  energy  lasers, 
the  window  heating  causes  a rapid  deterioration  of  ivs  transparency. 

In  fact,  at  present  there  are  no  existing  solid  windows  for  practical 
use  on  high  energy  lasers. 

A third  technique,  and  the  one  v/hich  Is  most  viable  for 
present  high  energy  lasers  of  all  types  — electrical,  chemical  and 
gasdynamic,  Is  tc  extract  the  beam  through  an  aerodynamic  window. 

To  understand  the  principle  of  such  windows,  recall  that  expansion 
and  shock  waves  In  supersonic  flows  can  support  substantial  pressure 
differences.  This  Is  Illustrated  In  Figure  81.  In  a GDL,  the  cavity 
pressure  Is  generally  0.1  atm,  and  the  outside  ambient  pressure  is  1 atm. 
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Referring  to  Figure  81,  expansion  and  compression  waves  can  clearly  act 
as  a boundary  between  such  pressure  ratios,  where  Pj  and  pg  represent 
the  cavity  and  ambient  pressures  respectively.  This  Is  the  underlying 
Idea  of  the  aerodynamic  windows  shown  In  Figure  82,  obtained  from 
Reference  79.  Here,  three  types  of  windows  are  Illustrated.  All  three 
use  a convergent-divergent  nozzle  to  expand  the  flow  to  supersonic 
speeds.  In  the  expansion  window,  an  expansion  wave  Is  formed  at  the 
nozzle  lip  on  the  cavity  side.  In  the  compression  window,  an  oblique 
shock  wave  Is  formed  at  the  nozzle  lip  on  the  ambient  side.  The 
shock -expansion  window  combines  both  of  the  above  wave  systems.  In  all 
three,  the  laser  beam  propagates  out  of  the  low  pressure  cavity, 
through  the  expansion  and/or  shock  wave zt  and  Into  the  ambient  air 
at  1 atmosphere.  Since  air  is  almost  transparent  at  10. 6p,  the  laser 
Intensity  is  not  attenuated.  Of  course,  with  the  aerodynamic  window, 
there  are  no  material  heating  problems. 

A drawback  of  aerodynamic  windows  is  that  they  require  high 
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mass  flow  rates  of  dry  air  or  Parmentier  has  shown  that  the  mass 

flow  Is  proportional  to  the  square  of  the  aperture  diameter,  D (see 

Figure  82).  His  results  are  shown  in  Figure  83  for  the  compression 

window,  and  in  Figure  84  for  the  expansion  and  shock-expansion  windows. 

Clearly,  for  meaningful  aperture  sizes  on  the  order  of  10-20  cm,  large 

79 

mass  flow  rates  are  required.  Parmentier  points  out  that  the  shock- 
expansion  window  requires  the  smallest  mass  flow  of  any  of  the  three 
types. 

There  are  density  variations  in  the  beam  direction  through 
an  aerodynamic  window,  hence  there  are  phase  distortions  introduced 
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Fig.  84  Mass  flow  required  for  expansion  window 
operation. 


by  the  window  Into  the  beam.  For  the  expansion  or  compression  windows, 

the  phase  distortion  Is  linear  which  results  In  simple  beam  shearing, 

but  no  degradation  In  the  far  field.  The  density  distributions  through 

a shock-expansion  window  calculated  by  the  method  of  characteristics 

79 

are  given  by  Parmentler  , and  are  shown  In  Figure  85.  The  corresponding 
phase  distortion  Is  shown  In  Figure  86. 

Parmentler  has  conducted  experiments  with  a shock -expansion 

t 

window  which  are  reported  In  Reference  79.  Among  these  tests,  an  Independent 
COg  laser  beam  from  an  electric  discharge  laser  was  passed  through  a 4 cm 
aperture  aerodynamic  window,  as  shown  In  Fig.  87.  The  results  for  the 
output  beam  Intensity  distribution  are  shown  In  rlgure  38,  which  Indicates 
only  a small  effect  of  the  flow  on  the  beam. 

i 

The  results  of  Parmentler  clearly  show  the  feasibility.  Indeed 
the  necessity,  for  aerodynamic  windows  on  high  power  gas  lasers. 
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IV  EPILOGUE 

It  Is  the  author's  hope  that  a serious  student  with  no  prior 
familiarity  with  gasdynamlc  lasers  can,  by  studying  these  Notes  and  by 
examining  the  cited  references,  bring  himself  to  the  forefront  of 
GDL  understanding  and  the  state-of-the-art.  However,  the  technology 
of  such  lasers  Is  developing  rapidly,  and  what  Is  the  state-of-the-art 
today  can  become  "classical"  tomorrow.  Therefore,  the  reader  Is 
encouraged  to  build  on  these  Notes  by  keeping  up  with  the  modern  literature 
and  Indeed  making  his  own  contributions  In  the  laboratory  or  In  the 
theory  of  lasers.  In  the  pressing  energy  conversion  problems  of  today, 
high  energy  lasers  will  have  a multitude  of  applications.  We  have 
just  scratched  the  surface  so  far. 
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ABSTRACT 

Recent  data  for  the  vibrational  energy  ex* 
change  rates  associated  with  COj-Nj-^O  gas  lasers 

are  compared  with  previous  compilations,  from 
this,  current  uncertainties  in  the  rate  data  are 
established.  Using  an  existing  computer  code,  a 
study  is  made  of  the  sensitivity  of  gasdynamic 
laser  performance  to  these  uncertainties.  The 
results  show  that  present*day  calculations  of 
gasdynamic  laser  gain  and  maximum  available  power 
are  subject  to  at  least  a 25  percent  Inaccuracy 
due  to  such  rate  data  uncertainties. 

I.  INTRODUCTION 

The  first  successful  demonstration  of  a 
CO2-N2  gasdynamic  laser  took  place  in  1965,  as 

described  by  Gerry^.  Since  then,  the  state-of- 
the-art  has  blossomed  and  matured  at  a rapid  rate, 
as  surveyed  in  References  2-4.  In  particular,  a 
number  of  gasdynamic  laser  experiments  have  been 

carried  out  in  arc  tunnel s^5*^,  shock  tunnel  s^'3^ 

and  combustion  driven  devices^*  In  many 

of  these  experimental  Investigations,  comparisons 
have  been  made  with  theory  for  small -signal  laser 
gain  and/or  laser  power  extraction.  Fair  agree- 
ment has  been  obtained  in  so.ite  but  not  all  cases. 
With  the  rapid  development  of  such  lasers,  there 
is  an  increasing  need  to  obtain  closer  agreement 
between  prediction  and  experiment;  indeed,  it  is 
desirable  to  calculate  gasdynamic  laser  perform- 
ance as  accurately  as  possible  — ultimately  to 
within  a few  percent. 

At  present,  there  Is  an  inherent  limitation 
to  the  accuracy  with  which  such  calculations  can 
be  made,  namely,  the  uncertainty  in  the  pertinent 
vibrational  kinetic  rates.  The  heart  of  a gas- 
dynamic laser  is  the  nonequilibrium  nozzle  expan- 
sion in  which  the  population  Inversion  is  created, 
and  the  laser  cavity  downstream  of  the  nozzle  exit 
wnere  power  is  extracted.  In  both  these  regions 
the  computed  nonequilibrium  flow  (and  hence  the 
theoretical  laser  performance}  is  dependent  on  a 
complex  finite-rate  vibrational  energy  exchange 
sechanism  which  is  still  not  yet  fully  understood. 
For  the  COg-^^O  mixture  common  to  most  gas- 
dynamic lasers,  at  least  eight  distinct  vibration- 
al kinetic  rates  must  be  utilized  in  the  calcula- 
tions. In  turn,  uncertainties  in  these  rates  will 
cause  Inaccuracies  in  the  computed  gasdynamic 
laser  performance.  Hence,  two  questions  Immedi- 
ately arise:  (1)  what  are  the  present  uncer- 


tainties in  the  measured  kinetic  rates,  and  (2) 
how  sensitive  are  calculations  of  gasdynamic  laser 
performance  to  these  uncertainties?  The  purpose 
of  the  present  investigation  is  to  answer  these 
questions. 

in  particular,  the  rates,  their  uncer- 
tainties, and  some  details  concerning  their 
effect  on  gasdynamic  laser  small -signal  gain  and 
maximum  available  power  are  discussed  in  Section  11. 
The  quantitative  sensitivity  of  gasdynamic  laser 
calculations  to  each  vibrational  kinetic  energy 
exchange  reaction  is  defined  in  Section  III,  and 
bar  charts  are  presented  which  graphically 
demonstrate  the  absolute  uncertainties  which 
currently  exist  in  such  calculations.  Finally, 
those  particular  reactions  which  have  the  most 
influence  on  gasdynamic  laser  gain  and  maximum 
available  power  are  clearly  delineated. 

The  calculations  of  gasdynamic  laser  perform- 
ance in  the  present  paper  were  made  with  an 
existing  computer  program  described  in  Reference 
13.  This  computer  code  is  based  on  the  time- 
dependent  nonequilibrium  nozzle  flow  analysis 
discussed  in  References  14  and  15,  and  applied  to 
gasdynamic  lasers  In  Reference  16. 

II.  RATES  AND  UNCERTAINTIES 

The  vibrational  kinetic  rates  pertinent  to 
C02*N2  gas  lasers  were  examined  and  compiled  by 

Taylor  and  Bittennan,  first  in  1967  (Reference  17) 
and  slightly  updated  in  1969  (Reference  lb). 

These  references  are  timely  and  excellent  surveys, 
and  have  become  a standard  source  for  laser  cal- 
culations. However,  in  1969,  References  17  and 
IB  clearly  show  tnat  a large  scatter  existed  in 
much  of  the  data;  moreover,  for  a few  rates,  only 
one  set  of  data  were  available.  Hence,  in  the 
present  paper,  additional  COj-^’^O  vibrational 

rates  measured  or  calculated  since  1969  are  pre- 
sented along  with  the  standard  Taylor  and  Bittermar 
data;  this  is  an  effort  to  assess  uncertainties 
which  currently  exist  in  these  rates.  (Some  of 
the  rates  measured  since  1969  have  been  recently 
catalogued  in  an  unpublished  memorandum  by  0.  W. 

Ilall^3^.  The  present  author  gratefully  ack- 
nowledges the  help  of  this  memorandum  as  a refer-  . 
ence  source.) 

The  reactions  that  must  be  Included  in  the 
calculation  of  small -signal  g\in  and  maximum 
available  power  in  gasdynamic  lasers  are: 


fThis  research  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under  AFOSR  Grant 
No.  74-2575,  with  Captain  Lloyd  R.  Lawrence,  Jr.  as  Program  Manager.  A portion  of  the  computer  time 
was  kindly  provided  by  the  Computer  Science  Center  of  the  University  of  (Maryland. 
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For  such  calculations,  the  pumping  rate  between  ex* 
cited  N2  and  C02(001),  as  well  as  the  Fermi 

resonance  between  C02  (100)  and  C02  (020),  are 

considered  to  be  sufficiently  rapid  that  local 
equilibrium  is  established  between  l<2  and  C02 

(v2)  as  well  as  between  C02  (v^)  and  C02  (v2). 

Hence,  the  rates  for  these  reactions  are  not  con- 
sidered here.  However,  It  should  be  noted  that 
recent  Information  has  raised  some  questions  about 
the  local  equilibrium  normally  assumed  between  the 
C02  (vj)  and  C02  (v^)  modes.  Specifically,  Rosser, 

Hoag  and  Gerry^2®^  present  data  which  indicate  that 
Such  equilibrium  may  not  exist  for  all  laser  con- 
ditions. This  contention  is  further  supported  by 

Bulthius^21^,  who  presents  a measured  rate  for 
energy  transfer  from  the  (100)  to  (020)  level 
which  is  an  order  of  magnitude  slower  than  previ- 
ously reported  values,  and  prompts  a corment  that 
the  vibrational  temperatures  of  the  C02  (vj)  and 

C02  (v2)  modes  are  not  always  the  same,  particu- 
larly during  peak  power  extraction.  On  the  other 

hand,  De  Temple  et  al^22^  countered  these  conten- 
tions by  presenting  measurements  for  the  decay 
rate  of  C02  (100)  which  are  an  order  of  magnitude 

faster  than  those  of  Rosser  et  al  ^20^.  De  Temple 
et  al  feel  that  their  measurements  are  more  repre- 
sentative of  the  Intrinsic  V-V  relaxation  process, 
whereas  those  of  Reference  20  are  more  character- 
istic of  a V-T  process  under  the  influence  of  a 
perturbing  radiation  field.  This  matter,  albeit 
an  Important  one,  is  still  in  a state  of  flux; 
until  it  is  resolved,  the  assumption  of  local 
equilibrium  between  the  C02  (vj)  and  C02  (v2)  modes 

appears  to  be  the  only  logical  recourse  for  gas- 
dynamic  laser  calculations. 

The  effect  of  HgO  on  the  laser  process  is 

governed  by  reactions  1-3  above.  The  role  of  H20 

is  quite  Important;  Indeed,  the  deactivation  of 
the  lower  laser  level  by  H20  has  been  considered  in 

the  past  to  be  the  pivotal  reaction  for  C02-N2-H20 

gasdynamic  lasers.  The  rate  data  for  this  reaction 
Is  compiled  in  Figure  1.  In  this  ?nd  other  figures 
to  follow,  the  lines  and  symbols  are  identified  as 
follows: 

(1)  The  vertical  bars  delineate  the  scatter 
of  experimental  data  shown  in  the  1969  paper  by 

Taylor  and  Bitterman^®* , 


(2)  The  solid  curves  delineate  the  rates  used 
by  the  present  author  in  earlier  calculations  of 

gasdynamic  laser  properties^2**’  23’  2*^. 

The  solid  curves  will  be  denoted  as  "reference* 
rates  In  the  present  paper.  Simple  correlations 
of  these  rates  are  given  in  Appendix  A of  Reference 
23,  which  is  a modified  version  of  the  earlier 
correlations  given  in  Appendix  A of  Reference  24. 
These  correlations  were  originally  obtained  from 
the  "best  fits"  given  in  the  1967  report  of  Taylor 

and  Bitterman'  which  for  some  reactions  differ 
from  the  results  shown  in  the  1969  paper  by  the 

same  authors ^ ®^. 

(3)  The  dotted  lines  represent  fast  and  slow 
rates  chosen  for  performing  numerical  uncertainty 
tests,  to  be  described  later.  In  those  figures 
where  only  one  dotted  line  appears,  the  reference 
curve  is  used  for  the  other  limit  on  the  rate. 

Note  that  the  dotted  curves  are  for  numerical 
experimentation  only;  they  do  not  necessarily 
reflect  the  actual  uncertainty  that  exists  in  the 
data. 


(4)  The  solid  and  open  symbols  are  new  data 
published  after  the  appearance  of  the  1969  paper 

by  Taylor  and  Bitterman^  . 

Finally,  a note  is  made  that  the  rate  con- 
stants and  their  effects  on  gasdynamic  laser 
performance  which  are  examined  in  the  present 
paper  are  also  discussed  in  depth  in  Reference  25; 
this  reference  snould  be  consulted  for  more  details. 

l.  co2  (v2)  ♦ h2d*±cc2  ♦ h2o 


The  rate  data  for  this  reaction  is  compiled 
in  Figure  1,  which  gives  the  constant  as  a function 
of  temperature.  The  wide  scatter  in  the  data 
demonstrates  the  large  uncertainty  which  still 
exists  in  this  rate;  indeed,  even  the  temperature 
dependence  of  the  rate  constant  is  not  clear.  The 

experimental  data  of  Bulthius  et  al^26,2^, 

Buchwald  and  Bauer^2®^,  and  the  calculations  of 

Sharma^2^  are  shown  along  with  the  original  spread 

of  data  from  Taylor  and  Bittorman^®^ . At  the 
present  time,  the  reference  rate  (solid  line) 
appears  to  be  as  good  a representation  of  this  rate 
as  any  other. 

Using  the  computer  code  described  in  Reference 
13,  gasdynamic  laser  calculations  have  been  made 
alternately  employing  the  fast  and  slow  rates 
shown  as  the  dotted  lines  in  Figure  1.  Some 
corresponding  results  for  small -signal  gain,  GQ, 

are  shown  as  the  dotted  curves  in  Figures  2 and  3. 

In  these  calculations,  only  the  rate  for  C02 

(v2)  ♦ h20?£C02  ♦ H20  is  varied;  all  other  rates 

are  held  fixed  at  their  reference  values.  The 
solid  curves  in  these  figures  are  the  results 
obtained  using  the  reference  rates  only.  In 
particular.  Figures  2 and  3 illustrate  the  peak 
small-signal  gain  as  a function  of  reservoir 
temperature,  T , for  the  following  cases  respective- 
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(1)  20:1  area  ratio  nozzle  with  throat  height 


h*  • 1.0  «m,  71  C02.  11  H20,  921  N2. 

(2)  50:1  area  ratio  nozzle  with  h*  ■ 0.365ron, 
71  COj.  3.51  H20  and  89.51  N2> 

Cases  (1)  and  (2)  are  typical  of  'first  and 
second  generation"  gasdynamlc  lasers  respectively, 
as  described  In  References  2*4.  In  Figure  2,  the 
results  for  the  fast  and  slow  rates  differ  apprec- 
iably. In  Figure  3,  which  is  a case  with  3.51  H^O 

content  in  the  gas  mixture,  the  difference  Is  not 
large  at  low  temperatures,  but  increases  at  higher 
temperatures.  For  this  case,  the  relatively  high 
water  content  promotes  rapid  equilibration  of  the 
lower  laser  level  with  the  translational  tempera- 
ture, and  tends  to  swamp  the  effect  of  a factor  of 
3 difference  between  the  fast  and  slow  rates. 
Additional  results  are  presented  and  discussed  in 
Reference  25. 

From  the  results  shown  In  rigures  2 and  3,  a 
Quantitative  measure  of  the  sensitivity  of  gas- 
dynamic  laser  gain  to  uncertainties  in  the  rate 
constant  can  be  established.  Such  quantitative 
"sensitivity  factors"  and  "absolute  uncertainties" 
are  discussed  in  Section  III  of  this  paper.  More- 
over, the  trends  shown  in  these  figures  clearly 

indicate  that  the  reaction  C02  (v2)  ♦ HgO  *2  C02 

♦ HgO  plays  an  important  role  in  the  gasdynamlc 

laser  process,  a fact  that  is  well  known.  However, 
there  are  other  reactions  which  are  equally  as 
important;  these  will  be  discussed  in  the  following 
paragraphs.  The  importance  of  all  the  reactions 
will  be  compared  in  Section  III. 

2.  n*  ♦ h2o  *2  :i2  ♦ h2o 


The  rate  constants  for  this  reaction  are  given 
in  Figure  4 as  a function  of  temperature.  The 
symbols  and  curves  have  the  same  meaning  as  before. 
The  post-1969  experimental  data  of  von  Rosenberg 

et  al^33^  and  Evans^'^  are  shown.  Comparing  all 
the  data  in  Figure  4,  a wide  discrepancy  exists 
at  the  higher  temperatures. 

The  effect  of  fast  and  slow  rates  on  peak  gain, 

gain  profiles,  and  maximum  available  energy,  e_  . 

max 

is  shown  in  detail  in  Reference  25,  and  will  be 
suimariaed  in  Section  III  of  this  paper.  However, 
from  these  results,  the  conclusion  is  made  that  the 

reaction  N2  ♦ H20  <2  N2  ♦ H.,0  is  very  important  to 

the  gasdynamlc  laser  process,  especially  for  the 
case  of  a 50:1  nottle  with  3.51  H20  (typical  of 

possible  second  generation  devices).  This  reaction 
is  someth!  ig  of  a “sleeper"  in  that  (to  the 
author's  krowledge)  its  Importance  has  not  been 
fully  appreciated  in  the  past. 

3.  C02  ( V3>  ♦ H20  S2C02**  {^)  ♦ HgO 


The  rate  constants  for  this  reaction  are 
plotted  in  Figure  5 as  a function  of  temperature. 

The  data  of  Rosser  and  Gerry^3^,  and  Heller  and 

Hoore^33^  are  compared  with  the  single  set  of  data 

from  Taylor  and  Bitternan^0^ . Again,  there  exists 
a definite  uncertainty  in  the  data,  even  with 


regard  to  the  correct  temperature  dependence. 

The  Influence  of  uncertainties  in  this  rate 
on  gasdynamlc  laser  gain  and  maximum  available 
energy  arc  snown  in  detail  in  Reference  25.  As 
In  the  case  of  the  previous  two  rates,  the  reaction 

• . **•  , 

C02  (vj)  ♦ HgO  *±C02  (v2)  ♦ H20  exerts  a strong 

influence  on  laser  performance,  as  will  be  dis- 
cussed in  Section  HI. 

4.  C02  (v2)  ♦ H2^tC02  ♦ N2 


The  rate  constants  for  this  reaction  are 
plotted  in  Figure  6.  For  this  rate,  Taylor  and 

Bitterman^*^  show  no  actual  data,  hence  no  bars 
are  shown  In  Figure  6.  Instead,  they  recommend  a 

rate  which  is  0.2  times  that  for  CO,  (v,)  * CO,  32 

fl3l  2 2 2 

2C02.  The  reference  rate1  ' , snown  as  the  solid 
line,  uses  a factor  of  0.25.  Honevar,  Oata  from 
Shields,  Warf  and  Bass^34^,  Merrill  and  Amme^35^, 

and  Bauer  and  Schotter^3®^  are  now  available,  and 
are  plotted  on  Figure  6.  From  this  recent  infor- 
mation, it  appears  that  the  rate  for  C02  (v2)  ♦ 
N23E±C02  ♦ «2  is  closer  to  0.5  times  that  for 
C02  U2)  ♦ C02«=^2C02.  As  will  be  snown  in  Sectir 

III,  this  rate  does  not  have  a strong  influence  on 
gasdynamic  laser  performance. 

5.  CO*  (v2)  ♦ C025±2C02 


The  rate  constants  for  this  reaction  are 
plotted  in  Figure  7 as  a function  of  temperature. 

The  recent  data  of  Eckstrom  and  Bershader^3^ 
are  also  shown.  This  rate  appears  to  be  well- 
known;  moreover,  as  will  be  shown  in  Section  III, 
it  has  a very  minor  Influence  on  gasdynamlc  laser 
performance.  Hence,  any  uncertainty  in  this  rate 
is  of  little  consequence. 

6.  CO*  (v3)  ♦ h2?±C02**  (v2)  ♦ N2 


The  rate  constants  for  this  reaction  are  given 
in  Figure  8 as  a function  of  temperature.  The  data 

of  Rosser,  Wood  and  Gerry^3®^  are  also  shown.  Then, 
appears  to  be  a reasonable  uncertainty  in  tnis  rate, 
which,  as  will  be  shown  in  Sectior.  Ill,  has  a strong 
influence  on  gasdynamic  laser  performance. 

7.  C02  (v3)  ♦ C025£  C0~*  (v2)  ♦ C02 


The  rate  constants  for  this  reaction  are  given 
In  Figure  9,  along  with  the  more  recent  data  of 

Reid  et  al^3^  and  Seery^.  The  uncertainty  is 
not  severe;  moreover,  as  will  be  shown  in  Section 
III,  the  laser  performance  is  insensitive  to  un- 
certainties in  this  rate.  Hence,  this  rate  is  of 
little  importance. 

8.  N*  ♦ I<2  32  2N2 


This  is  by  far  the  least  irportant  rate  in  the 
gasdynamic  laser  process,  as  will  be  shown  in 
Section  III.  Hence,  the  uncertainty  in  the  rate 
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shown  by  Tsylor  *nd  Bltterman*18)  Is  of  no  conse- 
quence in  laser  performance. 

III.  SENSITIVITY  OF  GASDYHAKIC  LASER 
PERFORMANCE  TO  UNCERTAINTIES  IN  RATES 

In  this  section,  the  aforementioned  rate 
proce:ses,  their  uncertainties,  and  their  effect 
on  the  gasdynamic  laser  process  are  quantified  and 
compared.  Specifically,  the  vertical  distance 
between  the  dotted  lines  (or  dotted  and  solid 
lines)  shown  in  the  previous  figures  for  the  rate 
constants,  k,  define  a "change  in  rate  constant", 
ak.  In  turn,  this  change,  Ak,  causes  a correspond- 
ing change  in  small-signal  gain,  aG0,  as  shown  for 

example  in  Figures  2 and  3.  Hence,  a sensitivity 
factor,  S,  can  be  defined  for  each  rate  as: 

S • (aG0/G0)/(Ak/k) 

Here,  ak  Is  specified  on  the  previous  figures,  and 
aCL  Is  calculated  from  the  results  of  gasdynamic 
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laser  calculations'  ',  as  shown  In  Figures  2 and  3. 
Therefore,  $ physically  represents  the  uncertainty 
In  gain  per  unit  uncertainty  in  rate  constant. 

The  sensitivity  of  gasdynamic  laser  gain  to 
the  various  rates  Is  clearly  shown  in  Figure  10 
for  a low  nozzle  area  ratio  case  typical  of  first 
generation  devices.  Here,  the  sensitivity  factor 
S is  plotted  as  a function  of  reservoir  temperature, 

T ; a curve  is  shown  for  each  of  the  rates  con- 
o 

sidered  in  the  present  investigation.  For  the  case 
shown,  the  laser  gain  is  most  sensitive  to  re- 
actions 1,  3 and  6,  and  least  sensitive  to  reaction 
8.  In  fact,  for  T • 1200"K  (which  is  typical  of 
first  generation  gasdynamic  lasers),  the  rather 
remarkable  observation  Is  made  that  the  reactions 

CO*  (v2)  ♦ H20  5±C02  + H20 

C02  (v3)  ♦ H20  5^  CO 2 (v2)  + rl20 

C02  (v3)  ♦ N23±co***  (v2)  ♦ :i, 

all  exert  equal  and  maximum  sensitivity  on  calcu- 
lations of  gasdynamic  laser  gain. 

Similar  results  are  shown  in  Figure  11  for  a 
case  typical  of  high  area  ratio,  hign  H.,0  content, 

second  generation  devices.  Here,  in  addition  to 
reactions  1,  3 and  G,  reaction  2 is  seen  to  be 
Important.  In  fact,  at  T0  ■ 1800°K  (which  is 

typical  of  second  generation  devices),  the  re- 
actions 

co2  (v2)  ♦ h2o  co2  ♦ h2o 
n2  ♦ h2o  u2  ♦ h2o 

C02  (v3)  ♦ H20  SiC02**  (v2)  ♦ h20 

exert  the  most  sensitivity  on  calculations  of  gas- 
dynamic laser  gain.  Reflecting  on  Figure  11, 
emphasis  is  made  of  the  Importance  of  the  reaction 

N2  ♦ H20  S?ii2  ♦ H20.  The  effect  of  this  reaction 

previously  has  not  been  fully  appreciated  by  the 
present  author,  nor  presumably  by  other  workers  in 
the  field  as  well. 


Considering  maximum  available  energy,  eWj(,  a 

sensitivity  factor  can  also  be  defined  analagous 
to  that  for  gain.  This  factor  is  shown  in  Figures 
12  and  13  for  the  first  and  second  generation  laser 
cases  respectively.  Curves  for  the  rates  involving 
deactivation  of  the  lower  laser  level  are  not  shown 
because  they  have  a minimal  effect  on  cmax* 

Comparing  Figures  12  and  13  with  10  and  11,  the 
reactions  that  exert  maximum  and  minimum 
sensitivities  on  e.,,.  are  the  same  that  e.:*rt 
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maximum  and  minimum  sensitivity  on  GQ. 

The  results  shown  in  Figures  10-13  are  un- 
certainties in  performance  per  unit  uncertainty 
in  rate  constant,  hence  they  Indicate  the 
sensitivity  of  gasdynamic  laser  calculations  to  the 
individual  rates.  However,  as  previously  discussed, 
some  rate  constants  are  known  more  accurately  than 
others;  hence,  the  question  is  now  posed,  what 

are  the  absolute  uncertainties  in  G„  and  e . due 

o max 

to  the  existing  absolute  uncertainties  in  rate 
constants?  These  existing  absolute  uncertainties 
are  given  below,  and  were  obtained  from  examin- 
ation of  the  data  scatter  shown  in  Figures  1,4-9. 

Reaction  Uncertainty  In 

data,  Ak/k 


1.  CoJ  (v2)  ♦ H203±C02  ♦ H20  3 

2.  fl2  ♦ H20£H2  ♦ H20  5 

3.  C02  (v3)  ♦ H20?C0***  (v2)  ♦ H?0  5 

4.  C02  (v2)  ♦ N2J?C02  + N2  4 

5.  C02  (v2)  ♦ C02^t2C02  0.5 

6.  C02  (v3)  ♦ N2?iC02  (v2)  ♦ N2  5.0 

7.  C02  (v3)  ♦ CO^tCO***  (v2)  ♦ C02  0.25 

8.  :i2  ♦ .'<2«t  2m2  0.5 


In  establishing  tne  values  of  tk/k  shown  above, 
some  weight  has  been  given  to  the  recent  data. 

In  addition,  in  reactions  2,  3 and  6,  which  effect 
the  deactivation  of  the  upper  laser  level,  more 
weight  has  been  given  to  the  data  scatter  at  the 
higher  temperatures  because  a large  part  of  the 
upper  level  deactivation  in  a gasdynamic  laser 
occurs  in  the  throat  region  of  the  nozzle,  where 
the  static  temperature  is  still  reasonably  high. 

When  the  sensitivity  factors  shown  in  Figures 
10-13  are  multiplied  by  the  above  absolute  rate 
data  uncertainties,  Ak/k,  the  overall  uncertain- 
ties in  and  e . are  obtained.  These  results 
0 max 

are  shown  as  bar  charts  in  Figures  14  and  15  for 
G0  and  emax  respectively.  In  these  figures, 

reservoir  temperatures  of  1200*K  and  1800*K  are 
chosen  as  typical  of  first  and  second  generation 
gasdynamic  laser,  respectively.  The  Important 
conclusion  apparent  from  Figures  TT  ana  "15  i~ 
that  gasdynamic  laser  calculations  of  gain'  and 
maximum  available  enemy,  based  on  present-day 
knowledge*  of~~ra tes  in  tne  existing- l iterature, 
ap£ca r to~be  suTTjcct  to"  at  least  a 25  percent 
uncertainty. 
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IV.  COMPARISON  WITH  EXPERIMENT 

Based  on  the  above  study,  the  rates  that  are 
used  In  the  existing  computer  code  of  Reference  13 
appear  to  be  reasonable,  with  the  exception  of 

Nj  ♦ H20  N2  ♦ H20  (Figure  4)  and  C02  (v2)  ♦ 

N«3£C02  ♦ N-0  (Figure  6).  In  Figures  4 and  6,  the 
dot-da$ned  curve  represents  a modification  to  the 
reference  rates;  these  modifications  are  prompted 
by  the  postal 969  rate  data.  When  these  two  rates 
are  Improved  as  shown,  the  computer  code  of 
Reference  13  yields  the  results  shown  In  Figure 
16,  which  Illustrates  GQ  as  a function  of  H20 

content  for  the  given  conditions.  Two  curves  are 
shown,  one  using  the  old  rates  and  the  other  using 
'the  improved  rates.  These  are  compared  with  ex- 
perimental data  obtained  from  the  shock  tunnel 
(41 1 

measurements  by  Vamos'  ' . The  good  agreement 
obtained  between  the  theory  utilizing  the  Improved 
rates  and  the  experimental  data  seems  to  sub- 
stantiate the  two  modifications  which  are  Independ- 
ently based  on  Figures  4 and  6. 

For  convenience,  the  rates  used  In  the  gas- 
dynamic  laser  computer  code  of  Reference  13, 
modified  by  the  two  Improved  rates  discussed  above, 
are  given  in  Appendix  A of  the  present  paper. 
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APPENDI X S 

The  vibrational  relaxation  times  used  for 
generating  the  theoretical  curve  (improved  rates) 
in  Figure  16  are  given  below.  They  are  the  same  as 
appear  in  Appendix  A of  Reference  23,  with  the 
exception  of  the  two  improvements  mentioned  in  the 
present  report.  The  nomenclature  is  the  same  as 
References  23  and  24. 

*Tap)c02-N2  - 1.3  x 105  (t*1/3)4*9 

(tap)C02-C02  * 0,27  (Tap)CQ2-N2 

^Tap^C02-h20  “ 5,5  x 10  2 
,09(tbp)N2-H2  " 93  (T'1/3)  * 4-61 

(Tbp);(2-co2  " *Tbp)/12-N2 

,09(tbp)N2-n2o  ’ 27-65  (T‘1/3)  ■ 3-2415 

109Mco2-co2  “ ,7-8  <t‘1/3>  - 1-808 

{tcP^C02-N2  " 2(tcp)C02-C02 


V.  CONCLUSIONS 

From  the  present  study  of  the  sensitivity  of 
gasdynamlc  laser  performance  to  uncertainties  In 
the  vibrational  rate  data,  the  following  conclu- 
sions are  made: 

(1)  Recent  data  have  been  compared  with  the 

rates  compiled  by  Taylor  and  Bittemian^8^  in  1969. 
This  comparison  shows  that  an  inordinately  large 
uncertainty  still  exists  In  many  of  the  important 
vibrational  rates  used  for  the  calculation  of 
C02  - i<2  - H20  gasdynamlc  laser  performance. 

(2)  As  a result  of  these  uncertainties, 
present-day  calculations  of  gasdynamlc  laser  gain 
and  maximum  available  energy  are  subject  to  at 
least  a 25  percent  inaccuracy. 

(3)  Because  the  existing  rate  data  are  not 
sufficient  to  confidently  calculate  gasdynamic 
laser  performance  to  within  a few  percent,  new 
efforts  should  be  made  to  more  accurately  measure 
the  CQ2-N2-H20  vibrational  exchange  rates.  Of 

particular  Importance  arc  the  rates  for 

co2  (v2)  ♦ h2o  x±co2  ♦ h2o 

U*  (Vj)  ♦ H20  5±C0~*  (v2)  ♦ H20 

co2  (v3)  ♦ n2  co2  (v2)  ♦ h2 

N*  ♦ h20  «£N2  ♦ I120 


log(rcp)m  u q « -20.4  (T*1/3)  ♦ 0.643  (For 

c ra2*H2°  TxGOO'-') 

logU.pU  u n ■ -20.4  (600'1/3)  ♦ 0.642 

c lu2-h2u  , constant  (For  MkTceOJV.) 

In  the  above  correlations,  (rp)  is  In  (w  see  - 
atm)  and  T is  in  "K. 
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fig ur*  2 - Diagram  of  • vupersonlo  airing  nozzle 
with  COj  injection  after  expansion  of  tj. 

The  aerodynamic  pioture  lmjrove*  u the  injection 
port  Is  Bond  upotrou*  I Mixing  is  Initiates  earlier, 
sad  the  shook  verse  lose  their  strength  | furtlior, 
the  loss  of  vibrational  snorter  in  the  mixing  cons 
should  act  Increase  excessively.  Hovrrer,  one  dose  not 
want  to  lnj  >at  C02  and  Be  far  upotreas  of  the  throat, 

sines  the  advantage  of  mixing  would  ho  lost.  Injection 
at  the  throat  or  slightly  tv  -natream  offers  a reasoma- 
bls  cooprccise  (Pig.  3).  In.ootlon  on  tha  aria  ia  pro- 
faned hero,  rather  than  injection  through  tho  valla, 
in  order  to  avoid  tho  presence  of  CO^  in  the  boundary 

layers  | tM«  is  of  Importance  if  tha  laser  boss  la  to 
tnvtree  these  boundary  layer,  aa  in  e multiple  para- 
llel nozzle  system  for  Instance.  A similar  system  has 
• (Q 

been  described  raoontly  by  Croshfco  end  covorkero4 


Figure  3 - Diagram  of  s nozzle  with  injection 
et  the  throat. 


Tha  injection  olot  is  located  et  e point  whore  the 
Koch  number  ie  of  the  order  of  1.?.  Cold  prorixud 

CC?  and  Ea  ere  Introduced  et  sonic  velocities  Into  tie 
h'2.  The  near  equality  of  the  preeouree  of  both  flow* 

vlll  minimize  tho  onrot  of  shock*  and  thoirlo*  wlMtlee 
permit  thnlx  mr.ooth  miring,  leva  to  tha  point  whore 
the  injcotlon  ie  made,  the  »2  coming  frou  tha  etsgna- 

tien  dumber  retains  about  of  lte  Initial  vilrn- 
tional  o nor  or  owing  to  its  lorg  relaxation  time  in 

tt*  eboonce  of  tha  Fa  and  C<>2  contaminants^10^.  In 

the  euboequoat  portion  of  the  flow,  the  largest  fr mo- 
tion of  the  eolllaional  deactivation  of  tho  h'2  la 

esuaed  by  the  C02  and  He,  in  the  immediate  vicinity 
of  the  lajeotor,  whore  the  atntio  temperature  la  still 


high.  A few  ceotiaeters  dovnetrees  of  tbs  throat, 
aost  of  tha  initial  >2  vibrational  saergy  baa  bees 

fro  mu  while  the  C02  + Be  end  the  have  eeeplsted 
their  airing.  The  1»2  has  shared  its  vibrational 
energy  with  sod*  Vj  of  C02,  causing  e high  popu- 
lation Inversion  and  crasting  a large  optical  gain. 

It  should  be  stressed  that  other  advantages 
are  alee  seen  la  the  Injection  of  room  temporeturs 
C02  end  Be.  firstly,  It  represents  a saving  la  initial 

enthalpy,  and  therefore  an  Increase  ia  therco- 
dynasio  efficiency  | aeeoadly,  aa  tbs  CO^  and  Ba 

expand,  their  temperature  is  lowered  further,  thus 
lowering  the  Kg  static  tempers ture  as  they  diffuse 

Into  it.  Smaller  area  ratio  noeiles  can  then  be 
selected  for  a given  final  atetlo  tempo -a ture. 

These  qualitative  understanding*  have  been 
substantiated  by  nv.no ri cal  calculations,  which  vers 
conducted  in  order  to  characterize  various  nozzle 
designs  in  tens  it  nixing  efficiency,  gain  cons- 
tant and  optical  power  avaibl*. 


ITI.  Purwrlcal  celculatloey 

Tbs  program  treats  tha  caa*  of  steady  state  tvo- 
dlDenalonal  turbulent  eupereonlo  idling  in  a fashion 

similar  to  the  Patankar  and  Spalding  method  for 
aubeonio  flows.  Kodlflcatlons  have  been  madein  order  to 
account  for  tha  euporsunlo  nature  of  tha  flove  and 
the  diverge non  of  the  nozzle,  hevsrthelase,  insta- 
bilities rsnetr  calculations  at  the  throat  iepos- 
tibia. 

Tne  equations  of  notion  include  conservation  of 
axial  momentum,  species,  end  total  enthalpy  | the 
terms  of  mon-'.-.tua,  species  and  enthalpy  transfer 
resulting  i rox  turbulence  are  ernroBsod  accordirg 

to  the  simple  Prandtl  theory  and  depend  on  ar. 
experimental  determination  of  a mixing  leng-.h  cons- 
tant. In  addition,  equations  are  provided  lor  ths 
vibrational  energy  exchange  bet  en  tha  various 
Boloculer  opecies.  For  that  purpose,  a ritraticnal 
model  inspired  from  the  model  proposed  by  H.C.  Ionov 

et  al  ^ for  the  CO,  and  K,  molecules  haa  been 

(ll)  2 2 

adopted  4 . Teres  vibrational  modes  ax*  considered, 

namely  V„  (for  Kg),  Vj  of  C02,  and  ths  group  of 

mods*  V]  and  'Jg  of  C02  which  are  known  to  t-e 

tightly  eouplod.  The  energy  exchange  botveer  tha 
vibrational  modes  and  translation  Is  described  by 
the  equations  of  r.ef.  13,  while  the  rate  oonatar.ts 
are  taken  from  Taylor  and  Eittorman'e  review 

paper  It  has  also  tees  assumed  for  sicpli-ity 

that  there  are  no  fluctuation*  in  these  terms  re- si- 
ting from  turbulcsco,  and  therefor*  the  mean  vain** 
of  tha  parameters  (concentration,  temperature,  etc.) 
can  be  taken  at  each  j articular  location. 

The  boundary  conditions  are  thooa  that  reign  tha 
euperr.or.lc  streams,  which  have  equal  static  pro  s-rts, 
and  vhdeb  come  into  contact  at  tho  injection  plane  i 
the  Koch  1 COg  ♦ Is  flow  1*  supposed  to  te  in  equi- 
librium at  room  temperature,  while  tha  roloe  1*7. 


C-2 


liuit},  static  end  vibrational  temperature*  for  the 
lj  flea  arc  supplied  bp  a progra*  developed  for 

yremdxed  1^  -CO^-  Be  flova  typioal  of  staadaid  CDL^U^ 

' mith  tbs  COj  sad  Be  cocoentraUana  being  set  at 

aero.  Integration  of  the  cos  Klees  transformed  partial 
dlffe rential  equations  is  then  carried  out  froa  the 
Injection  place  downstream,  with  an  implicit  nuaorlcal 
S chaos.  Tbs  eorputiag  ties  la  approximately  13 
■iau tea  on  ea  HK  360-50  computer  for  aa  axial 

(7) 

Slatanoc  of  15  oa.  Vote  that  Spenoor  and  ooeorlcore' 
■lap  corked  out  a similar  program  for  their  EP 
laser*  The  tain  differences  lie  la  their  erpliclt 
so base  and  their  constant  prosourt  assumption. 

Typical  results  for  tbs  arrangement  in  rig.  2 art 
presented  la  Pig.  4 and  5*  is  expected*  a large 
tup*  return  lucre  ass  la  found  in  the  mixing  region* 
and  the  gains  arc  small  and  unevenly  distributed. 


Plgu-  j 4 - Temperature  distribution  for  auparuonio 
cdxl  • of  a Knob  5 Ep  floe  at  500  4K  eith  a Kach 
1,  i > #K,  pure  CC.,  ^floe.  liota  the  x and  y scales 
are  difforent. 


• 20  «<cm~ 


Figure  5 - Cain  distribution  for  aujorornic 
mixing  of  a Fseh  5.75,  400°K  static  tor;«rntv.re 
and  19J0*K  vibratlcr.nl  teeporotura  floe  eith 

a 105»  COj,  SX#  Ee  floe  at  Koch  1 end  250°)'.,  Cain 

eae  enlculatod  for  the  P 30  line  (optimum  gain  io 
slightly  larger  for  lince  P 23  - P 2'). 

On  the  oontrr.ry,  large  gains  arc  eot  up  serose 
tbs  floe  for  ths  nrrangvnent  of  11  g.  3*  Ono  can  see 
(Tig.  6)  that  the  dintributlon  is  quite  uniform  and 
that  tbs  gain  dooe  not  fell  eff  along  ths  flee 
direotion*  Ths  calculatlorn  further  indicats  that  the 


free xing  efficiency  ^ approaches  604,  with  a mat 
efficiency  f of  d>< 


Plgure  6 - Cain  dintributlon  oa  the  P 16  line  for 
mixing  at  tba  throat.  Tba  stagnation  conditions 
for  II^  are  2500»I  and  15  atm.  Tba  30*  CO^,  70*  Be 

mixture  ls-at  room  teoperaturc  | the  half  slot 
height  of  the  injector  la  0.3  *su  The  throat 
height  for  the  *2  (between  the  Injector  and  one 

nozzle  wall  la  0.7  ea).  Floe  proportions  are 
44*  K2,  39*  He,  17*  C02« 

Hu  Prperlpental  : o suite 

Both  types  of  mixing  nozzles  were  tested,  i d.o. 

2 Ktf  arc  beater  provided  a 100  g/a  floe  of  Ji'2*  nt 

temperatures  ranging  froa  2000  to  40G0(Z,  Ve  failed 
to  dstoct  ary  appreciable  optical  gain  eith  the 
first  type  of  arrangement  In  a 4 cm  eido  tu-o— 
dimensional  ehannol,  as  had  been  anticipated. 

The  mixing  assembly  for  tla  otner  errsngeaont  is 
sketched  in  Fug.  7.  Tho  threat  l.slfht  for  the  K2  la 

approximately  O.f  ma  on  cither  aide  of  the  injector. 
Various  injectors  vare  tooted,  eith  O.B  cm  diameter 
injection  holes  orior.Ud  eitr.er  on  axis  as  in  Tig.  7 
(injector  1),  or  nltcrr.stoly  up  end  dovn  et  rn  angle 
of  22°  from  tr.e  xid  plane  for  a better  mlilrg 
(injector  II)  t all  the  holes  were  spaced  at  1.2  mm 
sport. 


figure  7 - Sketch  of  a Koch  5 norzle  with  C02 
Injection  at  tha  threat  j tho  Initial  nupon.onic 
exionnlon  ii.-le  1e  7C*.  Tr.e  nottle  hvi  a 10  cm 
long  diver, -v:.t  jroflle  jlua  another  10  cm  long 
parallel  uectlon. 
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As  mixing  of  tho  two  flow  mi  investigated. 
Samples  nn  taken  alsultanaously  it  various  positloua 
la  tbi  flow  and  ruoaequantly  analysed  rlth  ■ gas 
chromatograph.,  Figure  8 giv w • typical  p^oflli  versus 
position  y ecroaa  nets  la  for  lnjactor  II  it  a diatanoa 
X ■ <9  ea  lro*  tba  noaala  throat. 


y(nwn) 


Figure  B - CC>2  do  la  fraction  19  ca  downstream 

froa  tha  throat,  for  tba  Injection  at  the  throat 
(tba  Ea  concentration  la  proportional  to  that  of 
COj).  Tba  Kj  stagnation  teeperatura  la  2000"£, 

with  oean  proportions  of  4fi£  15^  CO^  and 

39*  Be. 

Ortdcnl  hor.-uvneitv 

Scbllarcn  pictures  of  tha  flow  have  been  token 
with  1 ^is  exposure  tines,  revealing  the  ebaonce  of 

•hocks  | the  turbulence  lavol  of  tha  nlxc-d  flow 
waa  not  found  tc  exceed  significantly  that  of  tha 
pur*  heated  h'2  with  tha  injector  renovod  free  tha 

throat. 

Cain  reaeuTvc nte 

Tha  gain  eeaeureicnta  were  carried  out  with  tha 
•qulpnont  described  in  Kef.  14.  Tha  roanurenonta 
war*  node  with  tha  probe  bona  parallel  to  the  nozzle 
throat,  et  x ■ 15  ca.  Ho  signif leant  changes  in  the 
profiles  could  bo  found  froa  z ■ 10  ea  to  x ■ 19  ca. 
The  lower  halves  of  the  pain  profiles  era  plottod  in 
ri«.  9 versus  the  distance  y froa  tba  nld-plnro.  Tor 
injector  I,  tha  experimental  points  (open  circles) 
compare  well  with  theory  (solid  curve),  assuming 
affective  throat  blights  of  0.8  ca  for  I?2  and  0,15  ca 

for  COj  ♦ Ea,  with  respective  stagnation  temperatures 

of  2000“E  ar.l  300*1.  The  nixing  constant  was  token 
aa  a fitting  painnnter,  and  a value  of  0.03  was 
aalcetnd  for  it.  Injected  constituent  proportions 
(H^/CO j/na)  are,  In  nolos,  (CC/l5/i‘5),  tut  local 

proportions  In  tha  flow  vary  froa  thene  figures.  A 
central  dip  appears  on  the  oxperiner.tal  gain 
distribution  (at  y « 0),  but  lets  rarfcodly  than  on 
tba  theoretical  profile  | thio  can  bo  attributed  to 
tha  dlffervir.ee  bo  twee  the  actual  injector  design, 
which  is  composed  of  a scries  of  holes,  and  tho  r intli 
•lit  injector  colei  taken  for  tha  calculation.  Tha 


actual  Injector  perwlta  Rj  to  paaatrata  Into  tba 
aid-plane  between  tba  injection  ho  las. 


Figure  9 - Cain  distributions  across  tha  noxtla, 
at  x » 15  cm  downstream  froa  tba 
throat. 

For  Injector  II  (heavy  dots),  one  can  see  that 
tba  mixing  la  bora  efflclent.Tba  data  (takas  for 
tba  coalitions  of  Figure  6)  conflrn  that  tha  cl  ring 
la  near  satisfactory  for  |y|<15  ca.  Tha  stagnation 
temperature  la  2000*]:,  tha  freezing  efficiency  la 
of  the  order  of  ec£,  mean  Injected  constituent 
proportions  (<6/l5/39)  being  adjusted  for  taxinua 
gain  | tho  stagnation  pres cures  are  9 end  14  ate  for 
Kj  and  C02  ♦ Ea  respectively.  Tha  slightly  higher 

stagnation  pressure  for  C02  ♦ Ea  perolta  deeper 
penetration  Into  Kp.  Bata  taken  with  tha  ease  Injector, 
but  with  a 3000*X  stagnation  temperature,  BO  g/s  R2 

Ease  flow  rate  and  s (45/12/43)  cirtura  are  presented 
In  Figure  1C.  Tha  large  gains  cne  observes  car.  be 
Balntriinrd  up  to  iOC\ : K trmpor&ture,  provided  the 
C02  Dole  fraction  la  lowered  to  ICb.  approximately. 


Figure  10  - Cain  distribution  for  a 3000*1  h'2 

stagnation  tr-p.crature  with  • lower 
C02  Dole  fraction. 

Oscillator  err^rtr.mta 

A laser  cavity  was  also  operated,  with  a stable 
resonator  configuration.  The  mirrors  had  to  be 
recerrod  eonavhat  from  the  flow,  in  order  to  avoid 
doiva  by  the  nur-’ruus  eolld  particult-a  froa  tha 


c-  < 


aro-b**t«r  carried  tor  B?.  Virloua  output  mirror*  nri  5. 

tasted,  including  G*  or  double  RaCl  flat*,  and  tool* 
coupled  reflector*.  A cv  pov*r  of  2 IV  v*x  obtained 
with  th*  C*,  In  • **11  collimated  2.5  ca  dla  beam, 
for  O a tarnation  temperature  of  2800#K.  Rough 
calculation*,  toeaad  on  th*  120  kV  enthalpy  t- d to  th*  “• 
>2  which  actually  flow*  through  th*  earl  toy  and 

eoDtrltout**  to  tha  optical  pov»r,  gin  an  affielsncy 
.of  1.6JC  and  a epoclflo  poor  of  25  kV  per  kg/a.  '.'he 9 a 
figure*  coepar*  **11  with  value*  of  0.4?  and  4 kV  par 

k&/art*p*ctlvel y obtained  toy  Cany 


V.  Conclusion 


Extremely  proclelng  results  hara  been  obtained 
dth  a C02  GDI.  Nevertheless,  a number  of 

axparlasnta  remain  to  be  don#  an  thl*  ays tern. 


1)  A combustion  powered  system  la  preferable,  elne* 

tha  an  heater  1*  not  very  efficient.  On*  might  g( 

M*  th*  ooabuatlon  of  end  B2  0^»  which  give*  t' 

a mixture  of  I2  CO  et  ' .00*1  for  2 1 1 

proportions.  Although  CO  not  erpected  to  ruin 

th*  efficiency,  eons  loos  cay  occur  through  t0a 

collisions!  deactivation  of  CO^  In  the  id. ring  cone. 

2)  Eran  though  the  He  mole  fraction  Is  quite  low,  it 
seems  prof arable  to  replace  it  by  H^O  or  H2  for 

reduced  operating  cost*. 

j)  The  stagnation  presnures  oust  be  rained  for 

diffuser  axiaust  Into  the  atmosphere.  With  a Mach 
number  of  5,  a otagnation  pressure  of  30  eta  le  jj, 

needed  for  K^. 

Finally,  this  rising  teciirdque  r.lpht  alno  prove 
useful  for  er.enlcal  mixing  and  electrical  nixing  45 

(7  8 16)  ** 

laser*  v ’ ' ',  by  allowing  larger  gas  pressures 

ind  better  optlonl  homogeneities. 
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